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ABSTRACT 

Motivated by the long-standing "luminosity problem" in low-mass star formation whereby protostars 
are underluminous compared to theoretical expectations, we identify 230 protostars in 18 molecular 
clouds observed by two Spitzer Space Telescope Legacy surveys of nearby star-forming regions. We 
compile complete spectral energy distributions, calculate L^oi for each source, and study the proto- 
stellar luminosity distribution. This distribution extends over three orders of magnitude, from 0.01 
Lq - 69 Lq, and has a mean and median of 4.3 Lq and 1.3 Lq, respectively. The distributions are 
very similar for Class and Class I sources except for an excess of low luminosity (-Lboi ^5 0.5 L©) 
Class I sources compared to Class 0. 100 out of the 230 protostars (43%) lack any available data in 
the far-infrared and submillimeter (70 /im < A < 850 /im) and have Lboi underestimated by factors 
of 2.5 on average, and up to factors of 8 — 10 in extreme cases. Correcting these underestimates for 
each source individually once additional data becomes available will likely increase both the mean and 
median of the sample by 35% - 40%. We discuss and compare our results to several recent theoretical 
studies of protostellar luminosities and show that our new results do not invalidate the conclusions 
of any of these studies. As these studies demonstrate that there is more than one plausible accretion 
scenario that can match observations, future attention is clearly needed. The better statistics provided 
by our increased dataset should aid such future work. 

Subject headings: stars: formation - stars: low-mass - stars: luminosity function, mass function - stars: 
protostars 



1. INTRODUCTION 

Low-mass stars form from the gravitational collapse of 
dense molecular cloud cores of gas and dust (e.g., Be- 
ichman et al. 1986; Di Francesco et al. 2007). During 
the collapse process material accretes from the core onto 
the protostar. In this paper the term protostar is used 

1 Department of Astronomy, Yale University, P.O. Box 208101, 
New Haven, CT 06520, USA 

2 michael.dunham@yale.edu 

3 National Optical Astronomy Observatories, Tucson, AZ, USA 

4 Department of Astronomy, The University of Texas at Austin, 
2515 Speedway, Stop C1400, Austin, TX 78712-1205, USA 

5 Herzberg Institute, National Research Council of Canada, 5071 
W. Saanich Road, Victoria, BC V9E 2E7, Canada 

6 Center for Interdisciplinary Exploration and Research in As- 
trophysics (CIERA), Department of Physics & Astronomy, 2145 
Sheridan Road, Evanston, IL 60208, USA 

7 Institute for Astronomy, University of Hawaii at Manoa, Hon- 
olulu, HI 96822, USA 

8 Department of Astronomy, University of Massachusetts, 
Amherst, MA, USA 

9 Astrophysics Group, Physics, University of Exeter, Exeter EX4 
4QL, UK 

10 Department of Astronomy, University of Maryland, College 
Park, MD 20742, USA 

School of Physics and Astronomy, Cardiff University, Queens 
Buildings, The Parade, Cardiff, CF24 3AA, UK 

12 Herschel Science Centre, ESAC-ESA, P.O. Box 78, 28691 Vil- 
lanueva de la Canada, Madrid, Spain 

13 HarvardSmithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138, USA 

14 Space Science Institute, 4750 Walnut Street, Suite 205, Boul- 
der CO 80301 

1 ° European Southern Observatory (ESO), Karl-Schwarzschild- 
Strasse 2, D-85748 Garching bei Miinchen, Germany 



to refer to the hydrostatic object at the center of a col- 
lapsing core. More evolved young stellar objects (YSOs) 
no longer embedded within and forming from their natal 
dense cores are not considered protostars by this defini- 
tion. 

Despite several decades of progress, many details re- 
lating to the accretion of material from dense cores onto 
protostars remain poorly understood. As mass accretes 
onto protostars the gravitational energy is liberated and 
radiated away as accretion luminosity. This luminosity, 
which depends on the mass accretion rate, current pro- 
tostellar mass, and current protostellar radius, can be 
used to study the mass accretion process and distinguish 
between different accretion models. 

Observational studies of protostellar luminosities are 
hindered by the fact that protostars are deeply embedded 
in dense cores, with most of their emitted luminosities re- 
processed to mid-infrared, far-infrared, and submillime- 
ter wavelengths by the dust in the cores. The first signifi- 
cant study of the protostellar luminosity distribution was 
presented in a series of papers by Kenyon et al. (1990, 
1994) and Kenyon & Hartmann (1995). They identi- 
fied 23 protostars in the Taurus- Auriga molecular cloud 
and calculated bolometric luminosities by integrating the 
observed spectral energy distributions (SEDs) using In- 
frared Astronomical Satellite (IRAS) 12 — 100 /im pho- 
tometry and longer-wavelength (sub)millimeter photom- 
etry from the ground, when available. They found that 
the protostellar luminosity distribution extended from 
0.09 - 22 Lq, with a mean and median of 2.3 L Q and 0.7 
Lq, respectively, and a strong peak around 0.3 Lq. 
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As first noted by Kenyon et al. (1990), their observed 
protostellar luminosities are lower than expected from 
simple theoretical predictions. Their argument is as fol- 
lows. First, they assumed that all observed luminosity is 
accretion luminosity, 



where f acc is an efficiency factor taken to be 1, M and R 
are the mass and radius of the protostar, and M acc is the 
accretion rate onto the protostar. By further assuming 
that the peak of the luminosity distribution is produced 
by low- mass stars with M = 0.1 M Q and R = 1 R©, 
they calculated an implied mass accretion rate of ~ 10~ 7 
Mq yr _1 . If some fraction of the observed luminosity 
arises from the protostar itself (contraction, deuterium 
burning, etc.), the implied mass accretion rate is even 
lower. 

In the simplest model, the collapse of a singular 
isothermal sphere initially at rest as first considered 
by Shu (1977) and later extended by Terebey, Shu, & 
Cassen (1984) to include rotation (often called the "stan- 
dard model" of star formation) , collapse proceeds in an 
"inside-out" fashion, beginning in the center of the core, 
moving outward at the sound speed, and giving rise to 
a constant mass accretion rate of M acc ~ 2 x 10 -6 M Q 
yr _1 for 10 K gas. This is over ten times higher than 
inferred by Kenyon et al. (1990), and will only scale up- 
ward as Macc oc T 3 / 2 for higher gas temperatures. Mod- 
ifications to the standard model, including non-zero ini- 
tial inward motions (Larson 1969; Penston 1969; Hunter 
1977; Fatuzzo, Adams, & Myers 2004), magnetic fields 
(Galli & Shu 1993a, 1993b; Li & Shu 1997; Basu 1997), 
and isothermal spheres that are not singular but feature 
flattened density profiles at small radii (Foster & Cheva- 
lier 1993; Henriksen, Andre, & Bontemps 1997) all tend 
to increase the accretion rate over that predicted by the 
standard model, making reconciliation between theory 
and the Kenyon et al. observations difficult. This has 
become known as the "luminosity problem." 

Identification of protostars and determining their lu- 
minosities were both greatly improved by the launch 
of the Spitzer Space Telescope (Werner et al. 2004) in 
2003. Many sites of star formation have been observed at 
wavelengths ranging from 3.6 to 160 /im through various 
Spitzer surveys. One such survey was the Legacy survey 
"From Molecular Cores to Planet Forming Disks" (here- 
after c2d; Evans et al. 2003), which observed 7 large, 
nearby molecular clouds and ^100 isolated dense cores 
and resulted in the discovery of very low luminosity ob- 
jects (VeLLOs), protostars with internal luminosities 16 
< 0.1 L embedded in dense cores (Young et al. 2004). 
Dunham et al. (2008) identified 15 VeLLOs in the c2d 
dataset, and detailed studies of several have confirmed 
their very low luminosities and status as embedded pro- 
tostars (Dunham et al. 2006; Bourke et al. 2006; Lee et 
al. 2009; Dunham et al. 2010b; Kauffmann et al. 2011). 

Both Enoch et al. (2009) and Evans et al. (2009) stud- 
ied the c2d protostellar luminosity distribution by using 
the Spitzer data to identify protostars, compiling com- 

16 The internal luminosity, L[ nt , is the luminosity of the central 
source and excludes luminosity arising from external heating. 



plete SEDs including far-infrared and (sub) millimeter 
photometry from the literature, and integrating these 
SEDs to determine Lboi- They found a total of 112 
protostars in the seven c2d clouds. Enoch et al. cal- 
culated mean and median values similar to those found 
by Kenyon & Hartmann (1995), but with their improved 
sample statistics they noted the presence of a larger frac- 
tion of sources at low luminosities (< 1.0 L©). Evans et 
al. (2009) included a correction for foreground extinction 
and calculated revised mean and median values of 5.3 L© 
and 1.5 L©. In a separate study, Kryukova et al. (2012) 
also derived the protostellar luminosity distribution for a 
number of star-forming clouds, including the c2d clouds. 
They found an even larger excess of low-luminosity proto- 
stars than found by Enoch et al. and Evans et al. Offncr 
& McKee (2011) argued that the higher observed lumi- 
nosities found when extinction corrections are applied, 
combined with a more realistic value of the efficiency fac- 
tor in Equation[T]of f acc ~ 0.5 to take into account both 
the powering of jets and winds and the effects of unseen, 
episodic accretion bursts, can essentially resolve the lu- 
minosity problem, although explaining the large fraction 
of sources at very low luminosities remains a challenge. 

Several recent theoretical studies have explored possi- 
ble resolutions to the luminosity problem, many of which 
were originally proposed by Kenyon et al. (1990). One 
possibility is that accretion is variable or episodic, with 
prolonged periods of low accretion punctuated by short 
bursts of rapid accretion. Numerous origins for such a 
process have been proposed, including gravitational in- 
stabilities in protostellar disks (e.g., Vorobyov & Basu 
2005, 2006, 2010; Machida et al. 2011; Cha & Nayak- 
shin 2011), a combination of gravitational and magneto- 
rotational instabilities in protostellar disks (e.g., Ar- 
mitage et al. 2001; Zhu et al. 2009a, 2009b, 2010), quasi- 
periodic magnetically driven outflows in the envelope 
(Tassis & Mouschovias 2005), decay and regrowth of 
MRI turbulence (Simon et al. 2011), close interaction 
in binary systems or in dense stellar clusters (Bonnell & 
Bastien 1992; Pfalzner et al. 2008), and disk-planet in- 
teractions (Lodato & Clarke 2004; Nayakshin & Lodato 
2011). Indeed, Dunham & Vorobyov (2012) showed that 
the iboi distribution predicted by the Vorobyov & Basu 
(2005, 2006, 2010) simulations, which feature highly vari- 
able accretion with episodic bursts, provides a reason- 
able match to the c2d observations presented by Evans 
et al. (2009). Alternatively, Offner & McKee (2011) pre- 
sented analytic derivations of the protostellar luminos- 
ity function for several different accretion scenarios and 
showed that accretion models that tend toward a con- 
stant accretion time rather than a constant accretion rate 
provide a good match to the Evans et al. c2d observa- 
tions. As a third alternative, Dalba & Stahler (2012) 
recently argued that external accretion onto collapsing 
cores from the surrounding background cloud will reduce 
accretion rates and luminosities. 

With 112 protostars spread over more than three or- 
ders of magnitude in ibob the c2d sample of protostellar 
luminosities is still somewhat limited by small number 
statistics. As a follow-up to c2d, the Spitzer Gould Belt 
Legacy Survey (hereafter GB; L. Allen et al. 2012, in 
preparation) observed most of the remaining clouds in 
the Gould Belt. In this paper we extend the identifica- 
tion of protostars and calculations of Lboi from Evans et 
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TABLE 1 

Molecular Clouds Surveyed by the c2d and GB Surveys 







Distance 


Distance 




Cloud 


Survey 


(PC) 


Reference a 


Data R,eference(s) b 


Aquila 


GB 


260 


Maury et al. (2011) 


Gutermuth et al. (2008); Maury et al. (2011) 


Auriga/California 


GB 


450 


Lada et al. (2009) 


H. Broekhoven-Fiene et al. (2012, in preparation) 


Cepheus 


GB 


200-325 c 


Kirk et al. (2009) 


Kirk et al.(2009) 


Chamaeleon I 


GB 


150 


Belloche et al. (2011a) 




Chamaeleon II 


c2d 


178 


Whittet et al. (1997) 


Young et al. (2005); Porras et al. (2007); Alcala et al. (2008) 


Chamaeleon III 


GB 


150 


Belloche et al. (2011a) 




Corona Australis 


GB 


130 


JNeuhauser & lorbrich (2008) 


Peterson et al. (2011) 


IC5146 


GB 


950 


Harvey et al. (2008) 


Harvey et al. (2008) 


Lupus I 


c2d 


150 


Comeron (2008) 


Chapman et al. (2007); Merin et al. (2008) 


Lupus III 


c2d 


200 


Comeron (2008) 


Chapman et al. (2007); Merin et al. (2008) 


Lupus IV 


c2d 


150 


Comeron (2008) 


Chapman et al. (2007); Merin et al. (2008) 


Lupus V 


GB 


150 


Comeron (2008) 


Spczzi et al. (2011) 


Lupus VI 


GB 


150 


Comeron (2008) 


Spezzi et al. (2011) 


Musca 


GB 


160 


Knude & Hog (1998) 


T. Huard et al. (2012, in preparation) 


Ophiuchus 


c2d 


125 


de Geus et al. (1989) 


Padgett et al. (2008) 


Ophiuchus North 


GB 


130 


Wilking et al. (2008) 


Hatchell et al. (2012) 


Perseus 


c2d 


250 


Enoch et al. (2006) 


j0rgensen et al. (2006); Rebull et al. (2007) 


Serpens 


c2d 


429 


Dzib et al. (2010, 2011) 


Harvey et al. (2006, 2007a, 2007b) 



a Reference for the distance quoted in this work. 

References presenting the Spitzer IRAC and MIPS observations. 
c Different regions within Cepheus arc located at different distances; sec Kirk et al. (2009) for details 



al. (2009) to the combined c2d+GB dataset. Our work is 
motivated by a desire for better underlying statistics in 
the observed protostellar luminosity distribution and im- 
proving the accuracy of the Lboi calculations by includ- 
ing additional data not yet available when the Evans et 
al. study was conducted. The organization of this paper 
is as follows: We describe our method in §21 including 
overviews of the c2d and GB surveys in §2.1[ the iden- 
tification of protostars in §2.21 the compilation of full 
source SEDs in H2.31 and the calculation of Lboi in ^2.41 
§3] summarizes our basic results. A discussion of these 
results is contained in §31 In particular, in §4j]we com- 
pare our results to the existing c2d f §4.1.1|) and Kryukova 
et al. (2012) (i j4.1.2jl results, in §4~2l we discuss several 
recent theoretical investigations of protostellar luminosi- 
ties, and in H4.3I we evaluate the accuracy of our L DO i 
measurements for sources with observed SEDs that are 
not well sampled in the far-infrared and submillimeter, 
and the effects of this incomplete sampling on our overall 
results. Finally, we outline important future work needed 
to further advance this topic in §3J and summarize our 
findings in §51 



2. METHOD 



2.1. 



Overview of the Surveys 

The Spitzer c2d survey (PI: N. J. Evans) conducted an 
imaging survey of seven large, nearby molecular clouds 
and about 100 isolated molecular cloud cores, and a spec- 
troscopic survey of selected targets. The science ques- 
tions motivating this survey and a summary of the obser- 
vation strategy are given by Evans et al. (2003). A sum- 
mary of the results from the survey of the large molec- 
ular clouds is given by Evans et al. (2009). The Spitzer 
GB survey (PI: L. E. Allen) was designed as a follow-up 
to the clouds portion of c2d and conducted an imaging 
survey of 11 nearby molecular clouds, completing most 
of the remaining clouds in the Gould Belt (L. Allen et 
al. 2012, in preparation; see also Gutermuth et al. 2008; 
Harvey et al. 2008; Kirk et al. 2009; Peterson et al. 2011; 
Spezzi et al. 2011; Hatchell et al. 2012). The two surveys 



obtained 3.6-8.0 fiin images with the Spitzer Infrared 
Array Camera (IRAC; Fazio et al. 2004) and 24-160 /jm 
images with the Multiband Imaging Photometer (MIPS: 
Rieke et al. 2004) of all 18 clouds. A standard pipeline 
developed by c2d was used for data reduction, source ex- 
traction, and band-merging to produce final source cat- 
alogs for both surveys and has been described in detail 
elsewhere (Harvey et al. 2006; Evans et al. 2007). 

Table [T] lists each cloud, the survey in which it was 
imaged (c2d or GB), the assumed distance to the cloud, 
the reference for the distance, and references of individ- 
ual studies of each cloud where the observation strategy 
and basic results are presented. These clouds were cho- 
sen to represent nearly all of the significant sites of star 
formation within the Gould Belt, with two major excep- 
tions: the Taurus and Orion molecular clouds. These 
two clouds were each the focus of separate, dedicated 
Spitzer Legacy surveys led by other groups, and fold- 
ing in their results with the c2d+GB clouds will be the 
focus of a future paper. The clouds listed in Table Q] 
span very large ranges of properties. For example, the 
total cloud masses range from a few hundred M Q (e.g., 
Chamaeleon II; Evans et al. 2009) to ~ 10 5 M (Au- 
riga/California Molecular Cloud; Lada et al. 2009), the 
star formation rates and star formation rate surface den- 
sities both span approximately two orders of magnitude 
(Evans et al. 2009; Heiderman et al. 2010), and the ra- 
tio of protostars to pre-main sequence stars (indicative 
of the amount of current star formation still on-going 
in the cloud) range from none (e.g., Lupus V and VI; 
Spezzi et al. 201l) to values in excess of 30% (e.g., Au- 
riga/California Molecular Cloud, Cepheus, IC5146, and 
Perseus; H. Broekhoven-Fiene et al. 2012, in preparation; 
Kirk et al. 2009; Harvey et al. 2008; Evans et al. 2009). 
We refer the reader to the individual cloud studies listed 
in Table [T] for more details and additional references. 

We caution that the distances to the 18 clouds sur- 
veyed are not all well-known, and some cloud distances 
are still under significant debate. One such example is 
the debate over the distance(s) to the Serpens and Aquila 
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regions. Recent VLBA parallax measurements led to a 
65% increase in the distance to Serpens compared to that 
assumed by the c2d team (429 vs. 260 pc; Straizys et 
al. 1996; Harvey et al. 2006; Dzib et al. 2010, 2011), and 
there remains debate whether or not Aquila is also lo- 
cated at this new, farther distance or even if all of Aquila 
is itself located at the same distance (e.g., Gutermuth et 
al. 2008; Maury et al. 2011). We do not list formal dis- 
tance uncertainties in Table [T] as such uncertainties are 
very poorly characterized in at least some clouds. Instead 
we refer to the references listed in Table Q] for detailed 
discussions on the various methods used to derive dis- 
tances and the uncertainties in these methods. Future 
distance revisions will require future revisions to the re- 
sults presented in this study. 

2.2. Sample Selection 

Our method for selecting protostars from the c2d and 
GB observations closely follows the selection method 
used by Evans et al. (2009) for the c2d clouds. We sum- 
marize the main points here and refer to Evans et al. for 
more details. 

The data reduction pipeline creates band-merged 
source catalogs incorporating 2MASS and Spitzer 1.25 - 
70 fim photometry for each cloud. Candidate young stel- 
lar objects (YSOs) are identified using a standard clas- 
sification method developed for the Spitzer c2d and GB 
projects. This method is described in detail in Harvey et 
al. (2007b) and Evans et al. (2007) and summarized in 
all of the publications presenting individual cloud studies 
listed in Table [TJ Briefly, this method uses the Spitzer 
SWIRE Legacy survey of the ELAIS Nl extragalactic 
field (Lonsdale et al. 2003), processed to simulate the 
sensitivity and extinction distribution of the clouds in 
the c2d and GB surveys, to determine the positions of 
galaxies in three different Spitzer color-magnitude dia- 
grams. Each source extracted in the c2d and GB cloud 
catalogs with infrared colors indicative of the presence of 
dust (sources with colors that can not be explained by 
extincted background stars) is then assigned an unnor- 
malized "probability" of being a galaxy or YSO based on 
its position in each color-magnitude diagram, its K— [4.5] 
color, whether it was found to be extended in either of the 
two shortest Spitzer IRAC bands (3.6 and 4.5 /^m), and 
its flux density at 24 and 70 /im. The color and magni- 
tude boundaries, along with the final boundary between 
candidate YSO and candidate galaxy in unnormalized 
"probability" , are set to provide a nearly complete elim- 
ination of SWIRE sources. We refer the reader to Har- 
vey et al. (2007b) for further details on this classification 
method. Similar classification methods have been pre- 
sented by other Spitzer studies of galactic star-forming 
regions (e.g., Gutermuth et al. 2009; Rebull et al. 2010; 
Kryukova et al. 2012). 

In total, we identified 3239 candidate YSOs in the 
18 c2d and GB catalogs. All sources were visually in- 
spected to remove residual contaminants, including re- 
solved galaxies misclassified as candidate YSOs and im- 
age artifacts identified as point-sources by the automated 
pipeline but lacking true point-source detections in one or 
more bands (see Evans et al. 2009 for details) . Follow-up 
optical spectroscopy of targets in Serpens presented by 
Oliveira et al. (2009) led to the identification and removal 
of 11 background giants with infrared excesses. We lack 



the data required to identify and remove such objects in 
the other clouds. Oliveira et al. (2009) found a contami- 
nation rate of 25% in their Serpens study. Serpens (and 
Aquila) are likely the worst cases due to their close prox- 
imity to the Galactic plane (spanning Galactic latitudes 
ranging from 2° to 10°), although Romero et al. (2012) 
recently suggested the contamination rate is at least as 
high in other clouds as well, and Hatchell et al. (2012) 
found that 27% of their sample of candidate YSOs in 
Ophiuchus North selected via the c2d criteria were likely 
to be background giants based on proper motion argu- 
ments. However, 80% of the contaminants identified by 
Oliveira et al. and 75% of the contaminants identified 
by Hatchell et al. are classified as Class III YSOs, thus 
even if the overall contamination rate is as high as 25% 
- 30%, our inability to remove these contaminants will 
not significantly affect this study since it is only focused 
on the subset of YSOs that are considered to be proto- 
stars. Finally, a few known YSOs missing from the list 
of candidate YSOs due to missing photometry at one or 
more Spitzer wavelengths caused by saturation or non- 
detections from being too deeply embedded were added 
by hand. 

The above process resulted in a final list of 2966 YSOs 
(since all 2966 sources passed visual inspection, we have 
followed the terminology used by Evans et al. [2009] and 
dropped the word "candidate" at this point). This is 
nearly a factor of three increase over the 1024 YSOs 
identified in the c2d clouds alone by Evans et al. (2009). 
Many of these YSO populations have already been pre- 
sented and discussed in detailed studies of individual 
clouds (see Table Q] for references) and in an analysis 
of the star formation rates and efficiencies of the c2d 
and GB clouds based on a preliminary version of the fi- 
nal YSO catalog (Heiderman et al. 2010). A complete 
analysis of the full YSO population, implications for star 
formation rates and efficiencies in the Gould Belt, and 
the evolution and lifetimes of YSOs will be presented in 
a forthcoming paper (L. Allen et al. 2012, in prepara- 
tion). Here we focus only on the observed luminosities 
of protostars. 

The final sample of protostars is identified from the list 
of 2966 YSOs by examining the full SEDs compiled for 
each source (see below) and selecting only those sources 
associated with at least one (sub)millimctcr detection at 
A > 350 /zm, resulting in a final sample of 230 protostars. 
This is identical to the procedure followed by Evans et 
al. (2009) except they used a cutoff wavelength of 850 
/im; we modified this to 350 /im because of the large 
increase in available data at this wavelength. No intrin- 
sic protostellar colors were assumed and no additional 
color criteria were imposed. This decision is motivated 
by numerous recent studies that have used dust radia- 
tive transfer models to show that protostars observed 
through outflow cavities can resemble more evolved Class 
II or Class III sources in the infrared (e.g., Whitney 
et al. 2003; Robitaille et al. 2006; Crapsi et al. 2008; 
Dunham et al. 2010a). By selecting all sources asso- 
ciated with (sub)millimeter detections we recover such 
sources and identify all YSOs that are associated with 
dense cores, although future follow-up observations are 
required to remove true Class II or III sources simply 
seen in projection against a dense core. 

By requiring a (sub) millimeter detection, our method 
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requires the availability of (sub)millimeter surveys cover- 
ing the full extents of the clouds surveyed by c2d and GB. 
This is not always the described in more detail 

in the next section below. The effects of this limitation 
will be discussed in detail in §4.1.21 where we compare 
to a recent study that used very different methods for 
selecting protostars and did not require (sub)millimeter 
detections. 

2.3. Constructing Full SEDs and Correcting for 
Extinction 

Similar to Evans et al. (2009), we compiled as com- 
plete SEDs as possible for each of the 2966 YSOs. In 
addition to the 2MASS and Spitzer 1.25-70 /im pho- 
tometry provided by the source catalogs, we included 
the following: (1) optical photometry, where available 
from the literature, (2) Wide-field Infrared Survey Ex- 
plorer (WISE) 12 and 22 /im photometry from the all- 
sky catalog 17 , (3) selected other ground-based optical 
and infrared data as compiled by the authors of the de- 
tailed studies of individual clouds (see references in Ta- 
ble [U, (4) Spitzer 160 /im photometry for sources de- 
tected and not located in saturated or confused regions, 
calculated using aperture photometry and aperture cor- 
rections as given by the MIPS Instrument Handbook 18 ; 

(5) SHARC-II 19 350 /xm photometry, when available, 
from a targeted survey of protostellar sources (Wu et 
al. 2007; M. M. Dunham et al. 2012, in preparation); 

(6) SCUBA 20 4 50 and 850 (im photometry, when avail- 
able, from the SCUBA Legacy Catalog (Di Francesco et 
al. 2008); and (7) other (sub)millimeter photometry from 
unbiased surveys of molecular clouds, where available. 

For the last item above, other (sub)millimeter pho- 
tometry from unbiased surveys of molecular clouds, we 
used photometry from the following surveys: (1) A 
MAMB02 21 1.2 mm survey of part of Aquila (Maury 
et al. 2011); (2) A LABOCA 22 870 /im survey of 
Chamaeleon I (Belloche et al. 2011a); (3) A SIMBA 23 
1.2 mm survey of Chamaeleon II (Young et al. 2005); (4) 
A LABOCA 870 /nm survey of Chamaeleon III (Belloche 
et al. 2011b); (5) A Bolocam 24 1.1 mm survey of Ophi- 
uchus (Young et al. 2006); (6) A Bolocam 1.1 mm survey 
of Perseus (Enoch et al. 2006); and (7) A Bolocam 1.1 
mm survey of Serpens (Enoch et al. 2007). 

Summarizing the above information, we have access to 
complete (sub)millimcter surveys for only 6 out of the 18 
clouds (Chamaeleon I, Chamaeleon II, Chamaeleon III, 

17 Available atjMr^T77Jrea|^acxaltechjedu/cgi-bin/Gator /npjwscaii 

18 Available at http:/ /irsa.ipac. caltech.edu/data/SPITZER/docs/ 
mips / mipsinstrumenthand book/ 

19 The Submillimeter High Angular Resolution Camera II 
(SHARC-II) is a 350 /an bolometer array operated at the Caltech 
Submillimeter Observatory (Dowell et al. 2003). 

20 The Submillimeter Common-User Bolometer Array (SCUBA) 
was a 450 and 850 fim bolometer array operated at the James Clerk 
Maxwell Telescope. 

21 The Max-Planck Millimeter Bolometer 2 (MAMB02) was a 
1.2 mm bolometer array operated at the IRAM 30-m telescope. 

22 The Large Apex Bolometer Camera (LABOCA) is an 870 ^tm 
bolometer array in operation at the Atacama Pathfinder Experi- 
ment telescope (Siringo et al. 2009). 

23 The SEST Imaging Bolometer Array was a 1.2 mm bolometer 
array in operation at the Swedish-ESO Submillimeter Telescope. 

24 Bolocam is a 1.1 and 2.1 mm bolometer array operated at the 
Caltech Submillimeter Observatory (Glenn et al. 1998) 



Ophiuchus, Perseus, and Serpens), plus a partial survey 
ofAquila and piecemeal coverage of other clouds from the 
SCUBA Legacy Catalog (Di Francesco et al. 2008). This 
incomplete (sub) millimeter coverage will affect both our 
luminosity calculations and ability to identify protostars, 
and these effects are discussed in detail in 34.31 and §5.11 
Finally, before using the SEDs to calculate bolometric 
luminosities, we correct the photometry for foreground 
extinction. We wish to only correct for the foreground 
cloud extinction and not the local extinction from the 
dense core itself, as in the latter case the extincted emis- 
sion is reprocessed to longer wavelengths and included in 
our observed SEDs. Determining the true line-of-sight 
extinction to a protostar from the foreground cloud is 
not a trivial task. Following Evans et al. (2009), we as- 
sign extinction values to all 2966 YSOs (a sample which 
includes the 230 protostars identified in this work) as 
follows: 

1. We adopt extinction values from the literature for 
Class II and III YSOs (classified via infrared spec- 
tral index; see Evans et al. 2009 for details) in- 
cluded in published optical studies. 

2. We de-redden the remaining Class II and III YSOs 
to the intrinsic near-infrared colors of an assumed 
spectral type of K7, found to be fairly representa- 
tive of the majority of Class II and III YSOs in 
the c2d clouds (Oliveira et al. 2009, 2010; see also 
Evans et al. 2009 for details). 

3. We de-redden all of the Class I and Flat spectrum 
YSOs (again classified via infrared spectral index) 
in each cloud using the mean extinction toward all 
Class II YSOs in that cloud. 

The extinction values adopted for each of our 230 proto- 
stars following this procedure are listed in Table [4] Most 
of the protostars in a given cloud have the same adopted 
extinction value since most protostars are classified as 
Class I or flat spectrum via their infrared spectral index, 
although some have different values since no intrinsic pro- 
tostellar colors were assumed by our selection criteria and 
thus some Class II YSOs are classified as protostars (see 
g221 above). 

Once the extinction values are assigned, we use these 
values combined with the Weingartner & Draine (2001) 
extinction law for Ry = 5.5 to correct the photometry for 
extinction. The choice of the Ry = 5.5 law rather than 
the Ry = 3 1 law is motivated by several studies showing 
?m ffiMF^fe^^n^^ dense regions 

in which stars form (e.g., Chapman et al. 2009). While 
we do caution that our approach is somewhat crude, it is 
the best that can currently be done and is significantly 
more reliable than ignoring the effects of extinction alto- 
gether. 

2.4. Calculation of Evolutionary Indicators 

Once we have constructed full SEDs as described above 
in 32.31 we use these SEDs to calculate the bolometric 
luminosities (ibol) and bolometric temperatures (Tboi). 
Lboi is calculated by integrating over all detections, 



^bol 



And 2 



S v dv 



(2) 



G 



The bolometric temperature is defined to be the temper- 
ature of a blackbody with the same flux-weighted mean 
frequency as the source (Myers & Ladd 1993). Following 
Myers & Ladd, Tboi is calculated as 



11 fn°° vS v dv 
T hol = 1.25 x ltr 11 J > , K 

Jo S » dv 



(3) 



Tboi can be thought of as a protostellar equivalent of 
T e ff for stars; Tboi starts at very low values (~ 10 K) 
for cold, starless cores and eventually increases to T e ff 
once the core and disk have fully dissipated. The inte- 
grals defined in Equations [5] and [3] are calculated using 
the trapezoid rule to integrate over the finitely sampled 
SEDs. To avoid model or fitting uncertainties and focus 
only on the observations themselves, we do not extrapo- 
late beyond the shortest and longest frequences at which 
data are available and we do not interpolate over missing 
data. Instead, we explore the effects of missing data on 
our Lboi calculations in M4.31 We calculate Lboi and Tboi 
twice, once with the original, observed photometry and 
once with the extinction-corrected photometry. 

3. RESULTS 

For each of the 230 protostars identified following the 
selection method described above, Table|4]lists a running 
index, the cloud in which the protostar is located, the 
Spitzer source name (which also gives the coordinates), 
the assumed Ay for extinction corrections, the infrared 
spectral index 25 (a), Tboi, and Lboi calculated from both 
the observed and extinction corrected photometry, and 
a flag indicating whether or not each protostar has any 
available data at 70 /xm < A < 850 /im (see §4.3[) . In Ta- 
ble H] the extinction corrected values are denoted as a', 
T bol , and L' hol to differentiate them from the observed 
values. Throughout the remainder of this paper we con- 
sider only the extinction corrected values and drop the 
primes for simplicity. We do not give uncertainties for 
the Lboi derived in this work. Statistical uncertainties 
calculated by propagating through the uncertainties in 
the observed fluxes are on the order of 10%, but the true 
uncertainties are dominated by incomplete sampling of 
the SEDs and are impossible to calculate for each source 
individually. These uncertainties will be discussed fur- 
ther in 34~3l 

Figure Q] shows the distribution of the extinction cor- 
rected values of Lboi for all 230 protostars in log space. 
With a minimum and maximum of 0.01 Lq and 69 L Q , 
respectively, this distribution extends over greater than 
three orders of magnitude. The mean and median are 
4.3 L Q and 1.3 L Q , respectively. These statistics are 
summarized in Table El Also listed in Table [2] are four 
dimensionless quantities calculated from the luminosity 
distribution: the standard deviation of log Lboi, the ratio 
of the median to mean Lboi, the ratio of the maximum 
to mean Lboi, and the fraction of protostars with Lboi 
< 0.1 Lq. These particular quantities are motivated by 
the recent theoretical study of protostellar luminosities 
by Offner & McKee (2011), to which we compare our 
results below in 

Figure [2] shows the Lboi distributions separately for 

25 The infrared spectral index, a, is calculated over all 2MASS 
and Spitzer detections from 2 — 24 fim (Evans et al. 2007). 




-10 12 

log Lbd W 



Fig. 1. — Histogram showing the distribution of extinction cor- 
rected Lboi for all 230 protostars in log space. The bins are 1/3 dex 
wide, and the error bars show the statistical (V~N) uncertainties. 
The solid vertical line shows the approximate completeness limit 
of 0.2 L for the c2d+GB sample. 
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Fig. 2. — Histogram showing Lboi distributions in log space with 
1/3 dex bins. The shaded histogram shows the distribution derived 
in this study for the 65 out of 230 objects in the combined c2d and 
GB samples classified as Class protostars. The dashed histogram 
shows the same thing, except for the 120 out of 230 objects clas- 
sified as Class I protostars. The classification is based on T bo i 
calculated according to Equation \3\ using the extinction-corrected 
photometry and the Class boundaries defined by Chen et al. (1995). 
The solid vertical line shows the approximate completeness limit 
of 0.2 L for the c2d+GB sample. 



Class and Class I sources. We have used Tboi, calcu- 
lated using Equation[3l to classify our sources, since Tboi 
is one of several commonly used indicators of class and 
evolutionary status (e.g., Dunham et al. 2008; Enoch et 
al. 2009; Evans et al. 2009; Maury et al. 2011). Follow- 
ing Chen et al. (1995), Class sources are selected with 
the criterion that Tboi < 70 K and Class I sources are 
selected with the criterion that 70 < Tboi < 650 K. In- 
spection of Figure [2] reveals that the peak and extent of 
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TABLE 2 
Luminosity Distribution 
Statistics 



Parameter Value 



Total Number 230 

Mean 4.3 a L Q 

Median 1.3 a L Q 

Minimum 0.01 Lq 

Maximum 69 Lq 

Standard Deviation of log 0.73 

Median / Mean 0.3 

Maximum / Mean 16.0 

Fraction < 0.1 L 0.07 



a As described in EI4.3I once far-infrared 
and submillimeter photometry becomes 
available for the 43% of the sample lack- 
ing any available data at 70 fim < A < 
850 /im, the mean and median will likely 
increase to approximately 5.8 and 1.8 
Lq, respectively. The effects of includ- 
ing such data on the overall distribution 
°f ^bol? an d thus on the other quantities 
listed in this Table, can only be investi- 
gated once such data arc available. 

the Lboi distributions are similar for Class and Class I 
sources. The distributions have mean (median) values of 
4.5 L and 3.8 L (1.4 L Q and 1.0 L©) for the Class 
and I sources, respectively. However, there is one signifi- 
cant difference in that there is an excess of low luminosity 
Class I sources compared to the Class population. For 
the Class I population, 36% have L DO i < 0.5 L Q , whereas 
for the Class population, only 20% have such lumi- 
nosities. A K-S test on the two distributions returns a 
value of only 0.04, demonstrating that the difference at 
low luminosities is statistically significant. These results 
are similar to those obtained by Enoch et al. (2009) for 
a smaller sample. Very recently, several extremely low 
luminosity, Class sources have been discovered in cores 
classified as starless based on Spitzer observations (e.g., 
Chen et al. 2010; Enoch et al. 2010; Pineda et al. 2011; 
Schnee et al. 2012; Chen et al. 2012), emphasizing that 
at least some of this difference may be due to a bias 
against the lowest luminosity Class sources in Spitzer- 
selected samples. This point is further emphasized by 
the fact that the excess of low-luminosity Class I sources 
occurs below our approximate completeness limit of 0.2 
L Q (see below), where any such comparisons are limited 
in utility. The true similarity of the Class and Class 
I Lboi distributions must be revisited once current and 
future surveys with Herschel and ALMA detect and char- 
acterize the full population of extremely low luminosity 
protostars. 

We emphasize that the results presented here are the 
observed bolometric luminosities of protostars, which are 
not the same as the intrinsic protostellar luminosities. 
Departures from spherical symmetry break the correla- 
tion between observed and intrinsic bolometric luminosi- 
ties, and external heating from the interstellar radiation 
field breaks the correlation between bolometric and in- 
ternal luminosity. Regarding the latter, external heating 
can add up to several tenths of a solar luminosity depend- 
ing on the local strength of the interstellar radiation field 
and the core mass available to be heated externally (e.g., 
Evans et al. 2001) and can dominate the observed Lbol 
for the lowest luminosity objects. A few specific exam- 



TABLE 3 



L bol AND L int 


FOR VeLLOS 




Source 


-kbol 


Z/jut Reference a 


L1014-IRS 


0.34 


0.09 


1 


IRAM04191-IRS 


0.13 


0.08 


2 


L1521F-IRS 


0.13 


0.05 


3 


L328-IRS 


0.18 


0.05 


i 


L673-7-IRS 


0.18 


0.04 


5 


L1148-IRS 


0.13 


0.08-0.13 


6, 7 



a References: (1) Young ct al. (2004); (2) Dunham 
ct al. (2006); (3) Bourke ct al. (2006): (4) Lee ct 
al. (2009); (5) Dunham ct al. (2010b); (6) Kauffmann 
ct al. (2005); (7) Kauffmann et al. (2011). 



pies of this point can be found in recent, detailed stud- 
ies of individual VeLLOs that use continuum radiative 
transfer models to separate internal and external heat- 
ing and determine the intrinsic Lj nt . The observed Lbol 
and model-derived Lj nt for six such sources are listed in 
Tabled For at least 5 and possibly all 6, the observed 
Lbol are above 0.1 L Q while the model-derived L mt are 
below 0.1 L Q , qualifying them as VeLLOs. As a con- 
sequence, the fraction of protostars with Lbol < 0.1 Lq 
reported in Tabled (0.07) does not imply that 7% of the 
sample are VeLLOs; many more VeLLOs with Lbol > 0.1 
Lq are likely present in the sample. 

We have decided not to attempt to correct our luminos- 
ity distribution for source inclination or external heating, 
since any such corrections would be model-dependent 
(and in the case of external heating would require de- 
tailed modeling of all low-luminosity protostars, a project 
far beyond the scope of this paper) . What we present are 
simply the observed bolometric luminosities (after cor- 
recting for extinction) . Theoretical studies that attempt 
to explain the observed protostellar luminosity distribu- 
tion must take these considerations into account. 

Since the relationship between the fluxes in the various 
Spitzer bands and Lbol depends not only on distance but 
also on the detailed spectral shape of each source, local 
strength of the external (interstellar) radiation field, and 
total core mass available to be heated externally, there 



is no one completeness limit for each cloud or for the 
full c2d+GB dataset. In a detailed search for and study 
of low luminosity protostars in the c2d survey, Dunham 
et al. (2008) found that the c2d data are sensitive to 
protostars with L int > 4 x 10~ 3 (<i/140pc) 2 L . With 
cloud distances ranging from 125 — 950 pc, the resulting 
luminosity sensitivites range from 0.003 — 0.18 L©, or 
0.003 - 0.04 L© if IC5146 is omitted. However, this sen- 
sitivity is for Li n t rather than Lboi; as discussed above, 
the two are not the same for low luminosity protostars, 
with Lboi equal to or greater than Li nt depending on 
the details of the external heating. In another study, 
Enoch et al. (2009) estimated completeness limits of -Lboi 
~ 0.01 — 0.05 L Q for protostars in the c2d clouds, al- 
though they emphasized that there was significant un- 
certainty in deriving such limits. Taking into account all 
of the above information, we conservatively estimate that 
our sample is only complete for Lboi > 0.2 L Q (the sensi- 
tivity limit for IC5146, the most distant cloud, using the 
Dunham et al. [2008] relation and assuming no external 
heating) , and mark this limit with a solid vertical line in 
all figures presenting histograms of Lboi- The existence 
of protostars below these limits will be discussed in £15.31 

4. DISCUSSION 

In this section we discuss our results in comparison 
to other observational and theoretical studies of proto- 
stellar luminosities. In £14.11 we compare to two recent 
determinations of the observed protostellar luminosity 
distribution, and in ^4.21 we discuss several recent theo- 
retical investigations of protostellar luminosities. Finally, 
in ij4.3l we discuss the effects of missing far- infrared and 
submillimeter photometry on our derived luminosities. 

4.1. Comparison to Other Observations 
4.1.1. Comparison to c2d Results 

Evans et al. (2009) identified 112 protostars in the 
c2d survey and calculated their observed bolometric lu- 
minosities. Our methods for identifying protostars, as- 
sembling complete SEDs, and calculating Lboi are very 
similar to theirs. All of their protostars are included 
in our study, but we have expanded to the rest of the 
star-forming clouds observed by the GB survey and thus 
increased the number of protostars from 112 to 230. 
We have also made three changes to the ancillary pho- 
tometry included when assembling complete SEDs: (1) 
we have included 12 and 22 /im photometry from the 
WISE all-sky survey, (2) we have included additional 
SHARC-II 350 fim photometry from a targeted survey of 
nearby, low-mass star forming regions (Wu et al. 2007; 
M. M. Dunham et al. 2012, in preparation) that was not 
yet available when Evans et al. (2009) completed their 
study; and (3) we have not included any IRAS photom- 
etry. The last change is motivated by the superiority of 
WISE 12 and 22 /im and Spitzer 70 /im data to IRAS 
12, 25, and 60 /im in essentially all cases, and the ex- 
treme confusion from both nearby sources and ambient 
cloud emission in the IRAS 100 /im data. 

Figure |3] compares the Lboi distributions from this 
work and from Evans et al. (2009). The new distribution 
obtained in this study has a similar shape and extent to 
the c2d-only distribution, except now with better statis- 
tics. The medians are also quite similar, with values of 
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Fig. 3. — Histogram showing Lboi distributions in log space with 
1/3 dex bins. The shaded histogram shows the distribution de- 
rived in this study for the 230 protostars in the combined c2d and 
GB samples (see Figure [T] for error bars). The dashed histogram 
shows the distribution for the 112 protostars in the c2d sample as 
derived by Evans et al. (2009). The solid vertical line shows the 
approximate completeness limit of 0.2 Lq for the c2d+GB sample. 
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Fig. 4.— Spectral energy distribution of NGC1333-IRAS2A, plot- 
ted as S v versus v in linear space. The light shaded area shows the 
result of integrating under the curve when the WISE 12 and 22 /im 
photometry is included, whereas the dark shaded area shows the 
extra amount amount added to the integral when no photometry 
is available between 8 and 70 /an, as was the case for Evans et 
al. (2009). 



1.3 Lq in this work and 1.5 Lq in the c2d-only sample 
(Evans et al. 2009). A K-S test on the two distributions 
returns a value of 0.33, indicating they are not signifi- 
cantly different. 

Despite their general similarities, the two distributions 
do have different means: 4.3 L Q in this study versus 
5.3 L in the c2d-only sample (Evans et al. 2009). The 
mean is strongly influenced by the highest luminosity 
sources, several of which were overestimated by Evans et 
al. (2009). To understand the cause of this overestimate, 
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we note that there are 14 sources in the Evans et al. sam- 
ple saturated at 24 //m with Spitzer and thus lacking any 
photometry between 8 and 70 /im. By including WISE 
12 and 22 /um photometry, which was not available to 
Evans et al., our updated sample fills in this gap. Figure 
H plots the SED for NGC1333-IRAS2A (source 144 in 
Table 0]), one of the 14 sources saturated in the Spitzer 
24 /im observations. The SED is plotted as S v versus v 
in linear space rather than the more typical vS v versus A 
in log space since the former is the space in which the in- 
tegral in Equation [2] is calculated. The light shaded area 
shows the result of the integral when the WISE 12 and 
22 /im photometry is included, whereas the dark shaded 
area shows the extra amount added to the integral when 
no photometry is available between 8 and 70 /jm. 

As clearly demonstrated by Figure 0] omitting pho- 
tometry between 8 and 70 /jm can lead to significant 
overestimates of Lboi- In the specific case of NGC1333- 
IRAS2A, Evans et al. (2009) measured L bo i = 76 L Q 
whereas we measure Lboi = 22 L Q with the WISE 12 and 
22 /im photometry included, a factor of 3.4 lower. Our 
measurement is consistent with previous measurements 
of Lboi for this source whereas the Evans et al. value 
is not (e.g., J0rgensen et al. 2002). For the 14 sources 
saturated at 24 pm, Evans et al. (2009) measured Lboi 
ranging from 2.6 to 76 L Q , with a mean and median of 
30 and 27 Lq, respectively. For those same 14 sources 
and with WISE 12 and 22 /im photometry included, we 
measure Lboi ranging from 1.4 to 63 Lq, with a mean 
and median of 20 and 13 Lq, respectively. Most of the 
decrease in the overall sample mean from 5.3 Lq to 4.3 
Lq is a result of correcting this overestimate for several 
relatively high luminosity sources. 

4.1.2. Comparison to Kryukova Results 

Recently, Kryukova et al. (2012) presented observed 
protostellar luminosity distributions assembled from 
Spitzer observations of 11 molecular clouds: the 7 c2d 
clouds (Chamaeleon II, Lupus I, Lupus III, Lupus IV, 
Ophiuchus, Perseus, and Serpens), Taurus, and 3 mas- 
sive star- forming clouds (Orion, Cep OB3, and Mon R2). 
In total they identified 727 protostars in these clouds. 
Figure [5] compares our results. 

The left panel of Figure [5] compares the Lboi distribu- 
tions from this work and from Kryukova et al. (2012). 
We use the contamination-subtracted Lboi distributions 
from Kryukova et al. for this comparison. The two distri- 
butions are generally quite similar, but since Kryukova 
et al. include three massive star-forming clouds in their 
sample, environmental effects may mask our ability to 
properly compare the two results. Thus, the right panel 
of Figure [5] compares the Lboi distributions from this 
work and from Kryukova et al., where now both samples 
are restricted to the clouds common to both samples (the 
c2d clouds). 

Inspection of the right panel of Figure [5] clearly shows 
that, for the same clouds, Kryukova et al. (2012) find an 
observed distribution of protostellar luminosities that is 
generally shifted to lower luminosities compared to our 
results, with a much lower mean (2.3 Lq versus 4.2 Lq 
in our sample) and much higher fraction of protostars 
with Lboi < 0.1 L Q (22% versus 7% in our sample). A 
K-S test on the two distributions returns a value of 0.01, 
verifying that the two distributions are statistically dif- 



ferent. Kryukova et al. do not compile complete SEDs 
to use in calculating bolometric luminosities. Instead, 
for all sources they identify as protostars, they calculate 
Lmir, the mid-infrared luminosity from their 2MASS and 
Spitzer 1.25 — 24 /im data, and a, the infrared spectral 
index calculated from 3.6 to 24 /im. For sources with a > 
0.3 common to both their sample of protostars and the 
Evans et al. (2009) sample, they then derive the following 
empirical relationship: 

^MI2: = (-0.466 ± 0.014 x log(a) + 0.337 ± 0.053) 2 

(4) 

where Lboi is from Evans et al. (2009). They use this 
relation to calculate Lboi for all protostars with a > 
0.3, and the value of this relation at a — 0.3 to cal- 
culate Lboi for all protostars with a < 0.3. At least 
some of the discrepancy between our results and those 
of Kryukova et al. may arise because we have made sev- 
eral changes to t he SED s used to calculate Lboi, as de- 
scribed above in 34.1.11 To examine this possibility, we 
re-derived the above empirical correlation using our new 
values of Lboi for the sources common to both our sam- 
ple and the Kryukova et al. sample, and obtained the 
following modification using a linear least-squares fit: 

= (_o.298 ± 0.046 x log (a) + 0.270 ± 0.013) 2 

Lboi 

(5) 

We illustrate the effects of this modification in Figure |H1 
which re-creates Figure 5 from Kryukova et al. (2012). 

Using our modified relationship between Lboi, Lmir, 
and a, the changes to the Kryukova et al. (2012) Lboi 
values range from an increase by a factor of 1.9 to a 
decrease by a factor of two, depending on a for each 
source. The mean change of all sources common to both 
samples is an increase by a factor of 1.5. Such a change 
can explain much of the difference in means in the two 
samples (4.2 Lq in this study compared to 2.3 Lq in 
Kyuokova et al.), but cannot fully explain the excess of 
low-luminosity sources. Instead, the remainder of the 
discrepancy between our results lies in source selection. 

Kryukova et al. identify 43 protostars in the c2d clouds 
not identified by us or by Evans et al. (2009). Figure [7] 
shows the Lboi distribution for these 43 sources. Most 
have Lboi < 1.0 Lq, and while the effects of these "extra" 
low-luminosity sources are significantly mitigated by sta- 
tistical contamination corrections included by Kryukova 
et al., their net effect is to cause an excess of low- 
luminosity sources compared to our results. The main 
difference between our method of selecting protostars 
and that of Kryukova et al. is our requirement of at 
least one detection at A > 350 /im to ensure associa- 
tion with dense cores. Since the 43 sources shown in 
Figure [7] are not in our sample, they are not associated 
with (sub)millimeter detections and thus not associated 
with known dense cores. By our definition of a protostar 
(see ijl]), these are not protostars. 

However, it is possible that at least some of these 
sources are in fact protostars associated with relatively 
low-mass cores not detected by the (sub) millimeter sur- 
veys we used to compile complete SEDs. This is sup- 
ported by the fact that many such surveys have rela- 
tively high completeness limits (for example, the 50% 



10 



0.20 - 



| | c2d+GB (this study) 
' i Kryukova et al. (2012) 



m 0.15- 

01 
D- 

■g 0.1 oh 



0.05 - 



0.00 



-1 1 

log L bol (Lg) 



0.20 



m 0.15 

Q_ 



o 0.10 

CO 



0.05 



0.00 



| | c2d (this study) 



1 c2d (Kryukova) 



-2-10 1 

log L b0] (y 



Fig. 5. — Histogram showing L^oi distributions in log space with 1/3 dex bins. In both panels, the solid vertical lines show the approximate 
completeness limits of 0.2 Lq for the c2d+GB sample. Left: The shaded histogram shows the distribution derived in this study for the 230 
protostars in the combined c2d and GB samples (see Figure [l] for error bars). The dashed histogram shows the contamination-subtracted 
distribution for the 727 protostars identified by Kryukova et al. (2012). Right: The shaded histogram shows the distribution from this work 
when only including the sources in the c2d clouds. The dashed histogram shows the contamination-subtracted distribution from Kryukova 
et al. (2012) when only including sources from the same c2d clouds. 




Fig. 6. — Re-creation of Figure 5 from Kryukova et al. (2012), showing «/ (iMlR/^bol) (left) an d ibol/^MlR (right) vs. log(o) for the 
protostars common to both our sample and the Kryukova et al. sample. Lbd is from this work whereas Lmir (the mid-infrared luminosity) 
and a (the infrared spectral s lope) are given by Kryukova et al. (2012). The solid lines show the best-fit relation from Kryukova et al., 
whereas the dotted line shows the modified relation derived using our new values of Lbol ( see t ex t for details). 



completeness limits for the Bolocam 1.1 mm surveys of 
Perseus, Serpens, and Ophiuchus are 0.8, 0.6, and 0.5 
M©, respectively; Enoch et al. [2008]). On the other 
hand, some sources may be contaminants masquerading 
in the sample. Kryukova et al. (2012) made a careful 
attempt to correct for such contamination in a statisti- 
cal sense, but the resulting corrections are highly uncer- 
tain and may have been underestimated. For example, 
they applied their protostar selection criteria to the cat- 
alog produced by the Spitzer SWIRE Legacy survey of 
the ELAIS Nl extragalactic field (Lonsdale et al. 2003) 



to estimate the contamination from galaxies and remove 
the effects of this contamination from their luminosity 
distribution. Such an estimate is a lower limit only be- 
cause it does not take into account the fact that galaxies 
in their cloud source catalogs are observed through the 
extra extinction of the cloud itself, reddening all galaxies 
and thus increasing the number of galaxies with their as- 
sumed colors of protostars. Furthermore, Heiderman et 
al. (2010) recently showed that many sources selected as 
Class I YSOs by their infrared colors lacked the presence 
of warm, dense gas and are thus not protostars (and may 
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Fig. 7. — Histogram showing the Lbol distribution in log space 
for the 43 protostars in the c2d clouds identified by Kryukova et 
al. (2012) but not by us or by Evans et al. (2009). The bins are 
1/3 dex wide. 

not even be YSOs at all; see Heiderman et al. for details). 
Since Kryukova et al. used similar color-based selection 
criteria, it is plausible that not all of their sources are 
actually protostars. Our (sub) millimeter detection re- 
quirement should remove such fake sources. 

Ultimately, we conclude that there are limitations to 
both our method and that used by Kryukova et al. 
We are limited by the availability and sensitivity of 
(sub)millimeter surveys whereas they are limited by un- 
certain statistical corrections for contamination. We pre- 
fer our method because it ensures a reliable sample, but 
the Kryukova et al. (2012) results emphasize that this re- 
liability may come at the expense of completeness. Some 
of the "extra" 43 sources in the c2d clouds they identify 
but we do not may in fact be real protostars. Quantify- 
ing how many is simply not possible until future surveys 
in the far-infrared and submillimeter become available, 
as discussed further in M 5 . 1 1 

4.2. Theoretical Investigations of Protostellar 
Luminosities 

With the new observed protostellar luminosity distri- 
butions derived from large Spitzer surveys (Evans et 
al. 2009; Kryukova et al. 2012; this work), several recent 
studies comparing the predictions of theoretical accretion 
processes to the observations have been published. In one 
such study, Dunham & Vorobyov (2012) coupled two- 
dimensional radiative transfer calculations with the nu- 
merical hydrodynamical simulations of Vorobyov & Basu 
(2005, 2006, 2010). These simulations predict accretion 
rates that both generally decline with time and feature 
short-term variability and episodic bursts caused by disk 
gravitational instability and fragmentation. Dunham & 
Vorobyov (2012) used the core, disk, and protostellar 
masses, radii, and mass accretion rates predicted by the 
simulations as inputs to their radiative transfer calcu- 
lations. They included the effects of external heating 
in their radiative transfer models, and calculated model 
SEDs at all inclinations from the beginning of collapse 
until the end of the embedded phase. They used these 
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Fig. 8. — Histogram showing Lbol distributions in log space with 
1/3 dex bins. The shaded histogram shows the distribution derived 
in this study for the 230 protostars in the combined c2d and GB 
samples (see Figure[T]for error bars). The dashed histogram shows 
the model distribution derived by Dunham & Vorobyov (2012) 
based on a set of hydrodynamical simulations predicting accre- 
tion rates that both generally decline with time and feature short- 
term variability and episodic bursts caused by disk gravitational 
instability and fragmentation. The solid vertical line shows the 
approximate completeness limit of 0.2 Lq for the c2d+GB sample. 

SEDs to calculate Lbol in the same manner as observers 
(Equation [2J at all inclinations and all timesteps. Fi- 
nally, they assembled a theoretical prediction of the ob- 
served luminosity distribution by calculating the fraction 
of total time the models spend at each Lbol, weighted by 
inclination and initial core mass. 

Dunham & Vorobyov compared their model luminosity 
distribution to the c2d observations presented by Evans 
et al. (2009) and showed that the models generally match 
the shape and spread of the observed luminosity distri- 
bution, indicating that the underlying variable accretion 
process predicted by the Vorobyov & Basu (2005, 2006, 
2010) simulations offers a possible resolution of the lumi- 
nosity problem. Figure |S] compares our new observations 
with the Dunham & Vorobyov (2012) models. As our re- 
sults are generally similar to Evans et al. (2009) in terms 
of the shape and extent of the observed luminosity distri- 
bution (see §4.1.1[) . except with improved statistics, the 
basic conclusions of Dunham & Vorobyov are unchanged. 
While the agreement is not perfect, the models provide 
a reasonable match to the observed luminosity distribu- 
tion, with a K-S test on the two distributions returning a 
value of 0.41. The only significant discrepancy between 
our observations and the Dunham & Vorobyov (2012) 
models is at Lbol ^ 0.2 Lq. As discussed in detail by 
Dunham & Vorobyov (2012), this Lbol regime is at or 
below the completeness limit of the c2d and GB surveys, 
rendering a proper comparison of observed and model 
Lbol impossible at such luminosities (see also §5.3j) . The 
variable accretion process predicted by the Vorobyov & 
Basu (2005, 2006, 2010) simulations and considered by 
Dunham & Vorobyov (2012) remains a valid solution to 
the luminosity problem. 

Offner & McKee (2011) have also presented a recent 
theoretical study of protostellar luminosities. They de- 
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rived analytic luminosity functions for several different 
accretion scenarios and compared their results to the c2d 
observations presented by Evans et al. (2009). In addi- 
tion to comparing to the mean and median of the ob- 
served distribution, they calculated four dimensionless 
quantities to characterize their luminosity functions and 
compared to the same quantities calculated from the ob- 
servations. Their study was our motivation for calculat- 
ing and tabulating the quantities listed in the last four 
rows of Table [2J although we emphasize that the last 
quantity presented in Table [21 the fraction of protostars 
with Lboi < 0.1 Lq, is not the same as the fraction of 
VeLLOs from Offner & McKee (2011). Their quantity 
was calculated as the ratio of sources with Lboi < 0.14 
L© to those with Lboi £1-4 L Q and was chosen to pro- 
vide a direct comparison to the observations presented 
by Dunham et al. (2008). 

Based on comparing their analytic luminosity functions 
to the c2d observations, Offner & McKee (2011) con- 
cluded that accretion scenarios that tend toward a con- 
stant accretion time rather than a constant accretion rate 
are better able to match the observed protostellar lumi- 
nosity distribution. Given the similarity of our observed 
luminosity distribution to that of the c2d sample, the 
new observations presented here do little to change the 
findings of Offner & McKee. Their conclusion is in gen- 
eral agreement with a series of investigations by Myers 
(2010, 2011, 2012), who derived analytic luminosity dis- 
tributions based on simple models of protostellar evolu- 
tion assuming constant protostellar birth rates, accretion 
from both the dense core and from the surrounding ambi- 
ent medium ( "core-clump" accretion, to use their termi- 
nology) , and accretion durations set by the assumption of 
an equally likely stopping time. They showed that such 
models, which predict accretion rates that increase with 
protostellar mass and thus tend toward a constant ac- 
cretion time rather than constant accretion rate, exhibit 
good agreement with observed protostellar luminosities. 

While the Dunham & Vorobyov (2012), Offner & Mc- 
Kee (2011), and Myers (2010, 2011, 2012) models have all 
succeeded in matching the observed distribution of pro- 
tostellar luminosities, they do so with models featuring 
very different accretion properties. Like the collapse of a 
singular isothermal sphere, which does not match obser- 
vations, the Dunham & Vorobyov (2012) models feature 
time-averaged accretion rates (averaged over the full du- 
ration of the embedded phase) that do not vary with the 
final mass of the protostar. In other words, higher mass 
protostars require more time to form than low-mass pro- 
tostars. Their solution to the luminosity problem is to 
invoke variaiblity and episodic accretion bursts as pre- 
dicted by the Vorobyov & Basu (2005, 2006, 2010) simu- 
lations, which reduce the accretion rates for most times 
and increase them during short-lived accretion bursts. 
On the other hand, Offner & McKee (2011) and My- 
ers (2010, 2011, 2012) solve the luminosity problem by 
adopting models with time-averaged accretion rates that 
increase with the final mass of a protostar, so that all pro- 
tostars form in about the same amount of time regardless 
of their final mass. These studies emphasize that there 
is more than one possible resolution to the luminosity 
problem, and future theoretical work is needed to better 
decipher the implications of protostellar luminosities for 
the underlying mass accretion process and to distinguish 



between these different accretion scenarios. In this work 
we have assembled a larger dataset to which such work 
should compare. 

One weakness of both the Offner & McKee (2011) 
and Myers (2010, 2011, 2012) models is that they do 
not predict observed protostellar luminosities. The ob- 
served luminosity of a protostar includes accretion lumi- 
nosity, photosphere luminosity, and external luminosity 
from heating by the interstellar radiation field, and is de- 
pendent on the physical structure and inclination of the 
source. Offner & McKee (2011) only include the two 
internal luminosity components (accretion and photo- 
sphere), and Myers (2010, 2011, 2012) include only accre- 
tion luminosity. Furthermore, neither set of models takes 
into account the effects of source structure and inclina- 
tion. The evolutionary radiative transfer models pre- 
sented by Dunham & Vorobyov (2012) do take all these 
effects into account and calculate observed luminosities 
before comparing to observations, but they can also be 
criticized for assuming a very simple physical structure, 
not allowing for any variation in this structure, adopting 
a fixed interstellar radiation field with no variation in its 
strength, attenuation, or spectral shape, and weighting 
by only one of several formulations of the stellar initial 
mass function. Future theoretical work must build on 
the foundations laid by these recent studies to properly 
compare theoretical predictions to observed protostellar 
luminosities. 

4.3. Effects of Including FIR/SMM Photometry 

The SEDs of embedded sources typically peak around 
100 - 300 /xm (e.g., Andre et al. 1999; Enoch et al. 2009; 
Dunham et al. 2008); accurate sampling in the far- 
infrared and submillimeter is thus necessary to ensure ac- 
curate measurements of Lboi- Missing photometry near 
the peaks of the SEDs will result in underestimates of 
Lboi- As described above in H2.31 we included Spitzer 
160 fim photometry for sources detected and not lo- 
cated in saturated or confused regions, SHARC-II 350 
fim photometry when available from a targeted survey 
of protostellar sources (Wu et al. 2007; M. M. Dunham 
et al. 2012, in preparation), and SCUBA 450 (im pho- 
tometry when available in the SCUBA Legacy Catalog 
(Di Francesco et al. 2008). Out of the 230 total proto- 
stars identified in this work, 130 (57%) include at least 
one photometry point at 160, 350, or 450 fim. The other 
100 (43%) lack any available photometry between 70 fim 
and at least 850 fim, and sometimes between 70 fim and 
1.1 mm. The last column of Table 0] indicates whether 
each protostar has any such photometry available. 

How much have we underestimated Lboi for the 43% 
of the sample lacking any far-infrared and submillime- 
ter photometry, and what effect does this have on the 
derived luminosity distribution? To address these ques- 
tions, the left panel of Figure [§] plots the Lboi distribu- 
tion separately for the sources with and without at least 
one observed photometry point at 160, 350, or 450 ^m. 
As expected, the Lboi distribution for sources without 
any far-infrared or submillimeter data is shifted to lower 
luminosities. To quantify the amount by which Lboi is 
underestimated for these sources, we took the 130 proto- 
stars with at least one such observed photometry point, 
removed all data between 70 and 850 fim, and recalcu- 
lated Lboi- The results are shown in the right panel of 
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Fig. 9. — Left: Lbol histograms for the 230 protostars from this study in log space with 1/3 dex bins. The solid, shaded histogram shows 
the Lbol distribution for the 130 protostars with at least one detection at 160, 350, or 450 fim, whereas the dashed histogram shows the Lbol 
distribution for the 100 protostars with no available photometry between 70 and 850 fim. The solid vertical line shows the approximate 
completeness limit of 0.2 Lq for the c2d+GB sample. Right: The ratio of Lbol calculated with all available photometry to Lbol calculated 
with all detections between 70 and 850 fim removed for the 130 protostars with at least one detection in this wavelength range, plotted 
versus Tb i. 



Figure HI which plots the ratio of Lbol calculated with all 
available photometry to L DO i calculated with all data be- 
tween 70 and 850 /im removed versus Tboi- The mean and 
median of this ratio are 2.6 and 1.5, with a few sources 
showing underestimates in Lbol up to factors of 8 — 10. 
There may be some evidence for a trend of larger Lbol 
underestimates for more deeply embedded sources (those 
with lower Tboi). Such a trend is not surprising since a 
greater fraction of the total luminosity is emitted in the 
far-infrared and submillimeter for more deeply embed- 
ded sources, but, if present, the trend is not very strong. 
Even sources with a Tboi of several hundred K can have 
Lbol underestimated by a factor of 2 or more. 

Another way of examining the effects of missing far- 
infrared and submillimeter data is to use Herschel 100 — 
500 /im photometry. The Herschel Gould Belt survey 
is currently surveying all of the clouds included here 
and will eventually provide such photometry (Andre et 
al. 2010). Maury et al. (2011) recently used data from 
this survey to study the protostellar population of part of 
Aquila. The ideal comparison would be to include Her- 
schel photometry in our SEDs and recalculate Lbol for 
each source, but we are unable to do this since Maury 
et al. did not publish the Herschel fluxes. They did, 
however, publish their own measurement of Lbol for each 
source including the Herschel data. We have identified 
17 protostars in Aquila that overlap with their sample; 
in all 17 cases we lacked any far- infrared or submillime- 
ter photometry. Figure [TU] plots the ratio of Lbol from 
Maury et al. (2011) with the Herschel data included to 
Lbol from this work without the Herschel data versus 
Tboi from Maury et al. (2011) with the Herschel data 
included. The results are strikingly similar to our above 
analysis: including the Herschel far-infrared and submil- 
limeter photometry increases Lbol by a mean and median 
of 2.6 and 1.6, respectively, and by up to factors of 8 — 10 
in some cases. Furthermore, there is again some evidence 
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Fig. 10. — The ratio of Lbol with Herschel photometry to Lbol 
without Herschel photometry versus Tb Q i calculated with Herschel 
photometry for the 17 protostars in Aquila. Lbol and Tb D i with 
Herschel photometry are taken from Maury et al. (2011) and Lbol 
without Herschel photometry is calculated in this work with the 
observed photometry (see text for details). 

of a trend with Tboi, but the sample size is too small for 
robust statistics. This analysis should be repeated for 
a larger sample and with Lbol calculated in a consistent 
manner with and without the Herschel photometry, but 
the agreement with the above analysis is encouraging. 

Based on two separate investigations, we conclude that 
lacking any far-infrared or submillimeter photometry will 
lead to underestimates in Lbol by a factor of about 2.5 
on average, but up to factors of 8 — 10 in extreme cases. 
These findings are consistent with earlier results by Dun- 
ham et al. (2008) and Enoch et al. (2009), who also inves- 
tigated this topic. Since nearly half (43%) of our sample 



14 



currently lacks any such photometry, including Herschel 
photometry when it becomes available will affect our de- 
rived luminosity distribution. To determine how signifi- 
cant this effect will be, we increased the calculated Lboi 
for all 100 sources lacking any far-infrared or submillime- 
ter photometry by factors of 2.5 and found that the mean 
and median of the full sample of 230 sources increase to 
5.8 and 1.8 L©, respectively. Compared to the values 
listed in Tabled both quantities increase by 35% - 40%. 
Testing this conclusion and determining the full effects 
on the shape and extent of the luminosity distribution 
can only be done once the Herschel Gould Belt survey is 
complete and the photometry is available. 

5. FUTURE WORK 
5.1. Completeness and Reliability of the Sample 

In order to be classified as a protostar, we require at 
least one detection at A > 350 /im to signify the pres- 
ence of a dense core. However, as described in ^2.31 
we only have access to complete (sub) millimeter surveys 
of 6 out of the 18 clouds (Chamaeleon I, Chamaeleon 
II, Chamaeleon III, Ophiuchus, Perseus, and Serpens), 
plus a partial survey of Aquila (Maury et al. 2011) and 
piecemeal coverage of other clouds from the SCUBA 
Legacy Catalog (Di Francesco et al. 2008). We are al- 
most certainly missing protostars from the clouds where 
we have incomplete (sub)millimeter coverage. Further- 
more, where we do have coverage, the surveys are typ- 
ically incomplete to low-mass cores with M < 0.5 M© 
(e.g., Enoch et al. 2008). We are also likely missing some 
protostars in low-mass cores. 

Kryukova et al. (2012) have also studied the proto- 
stellar luminosity distribution. They do not require a 
(sub)millimeter detection to be included in their sam- 
ple. In a detailed comparison between our results and 
their results (see H4.1.2p we showed that they identify 
approximately 40% more protostars than we do in the 
same clouds and that these "extra" protostars mostly 
have Lboi < 1.0 L©. Some may be true protostars lo- 
cated in low-mass cores below the sensitivities of the 
(sub)millimeter surveys available to us, but some may 
also be contaminants masquerading as protostars in their 
sample. 

Based on the presently available data, we are unable to 
evaluate how significantly our sample lacks completeness 
and how significantly the Kryukova et al. (2012) sample 
lacks reliability. Both should be fully resolved once the 
JCMT SCUBA-2 Legacy (Ward-Thompson et al. 2007) 
and Herschel (Andre et al. 2010) Gould Belt Surveys are 
complete, since the two surveys together will fully char- 
acterize the population of dense cores in all of the Gould 
Belt clouds with sensitivities below 0.1 M©. Future work 
must revisit the topic of defining the sample of protostars 
once the results from these surveys are available. 

5.2. Accuracy of Lboi Calculation 

As discussed in detail in H4.31 nearly half (43%) of our 
sample lacks any detections between 70 /jm and at least 
850 /xm, and sometimes between 70 ^m and 1.1 mm. We 
have shown that such sources will have their Lboi under- 
estimated by a factor of about 2.5, on average, and up 
to factors of 8 — 10 in extreme cases. We have estimated 
that correcting these underestimates will increase both 



the mean and median protostellar luminosity by about 
35% - 40%. Verifying this estimate and examining the 
full effects of these underestimates on the shape and ex- 
tent of the luminosity distribution must be revisited once 
the Herschel Gould Belt survey is complete and provides 
100 — 500 /im photometry for our full sample. 

5.3. Extending to Lower Luminosities 

As discussed in fj3j our protostellar sample is generally 
only complete for Lboi > 0.2 L©, with the exact com- 
pleteness limit different for each source depending on its 
distance, detailed spectral shape, and amount of exter- 
nal heating. Does there exist a population of protostars 
with luminosities below the completeness of the c2d and 
Gould Belt surveys? 

Recent evidence suggests that the answer is yes. First, 
in a detailed study of the population of low luminos- 
ity protostars in the c2d clouds, Dunham et al. (2008) 
noted that protostars were found all the way down to 
the sensitivity limit and suggested that the lower limit 
to protostellar luminosities had not yet been found. Sec- 
ond, very recent work has identified a population of 
cores originally believed to be starless and undetected 
in c2d and other Spitzer observations of similar depth 
but found to be driving molecular outflows through sen- 
sitive (sub)millimeter interferometer observations. To 
date, seven such objects have been identified, and all but 
one are located in the Perseus Molecular Cloud (Chen et 
al. 2010; Enoch et al. 2010; Dunham et al. 2011; Pineda 
et al. 2011; Schnee et al. 2012; Chen et al. 2012; Pezzuto 
et al. 2012). 

Both the true number and evolutionary status of these 
objects remains unknown. Most have been discovered 
serendipitously through (sub) millimeter interferometer 
detections of outflows in observations targeting cores be- 
lieved to be starless. A full survey for such objects has 
not been possible to date due to the prohibitively large 
time requests that would be required, although such sur- 
veys should be possible in the very near future with 
ALMA. None of the sources were detected in the infrared 
with c2d or other Spitzer surveys of similar depth, im- 
plying upper limits of 0.01 - 0.1 L© for both L mt and 
Lboi- One source was detected in very deep, targeted 
Spitzer 70 /jm observations and found to have Lj nt ~ 0.01 
L© and L bo i ~ 0.2 L© (Enoch et al. 2010). They have 
all been proposed as candidate first hydrostatic cores, 
a short-lived stage between the starless and protostellar 
phases (Larson 1969). First cores can drive outflows, al- 
though there is still debate about the physical properties 
of such outflows and which, if any, of the sources have 
outflow properties consistent with theoretical predictions 
(e.g., Machida et al. 2008; Dunham et al. 2011; Price et 
al. 2012). Future work must be devoted to determining 
how many of these objects exist, their true evolutionary 
status, and the effects their existence will have on the 
protostellar luminosity distribution. 

6. SUMMARY 

In this paper we have studied the protostellar luminos- 
ity distribution based on data from two Spitzer Legacy 
surveys of nearby star-forming regions. We summarize 
our main results as follows: 

1. Starting from a list of 2966 Young Stellar Objects 
identified via their positions in various color-color 
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and color-magnitude diagrams, we identify 230 pro- 
tostars in the Spitzer c2d and Gould Belt Legacy 
surveys based on association with a dense core de- 
tected at (sub)millimeter wavelengths. We compile 
as complete SEDs as possible for all 230 sources, 
and use these SEDs to calculate Lboi and Tboi- 

2. The protostellar luminosity distribution extends 
over three orders of magnitude, from 0.01 L Q - 69 
Lq, and has a mean and median of 4.3 L Q and 1.3 
Lq, respectively. Several dimensionless quantities 
characterizing the shape of the distribution are also 
calculated and tabulated. 

3. The luminosity distributions are generally similar 
for Class and Class I sources, with mean (me- 
dian) values of 4.5 L Q and 3.8 L Q (1.4 L© and 1.0 
L Q ) for the Class and I sources, respectively. The 
only difference is an excess of low luminosity Class 
I sources compared to the Class population: 36% 
of the Class I sources have Lboi < 0.5 L compared 
to only 20% for the Class population. A K-S test 
confirms that this difference is statistically signifi- 
cant. 

4. Our derived luminosity distribution is similar to 
that obtained by Evans et al. (2009) from the c2d 
data, except with better statistics. The most sig- 
nificant change is that we have added additional 
data to improve the accuracy of the Lboi measure- 
ment for sources saturated at 24 /jm with Spitzer, 
reducing Lboi by factors of ~ 2 — 3 for 14 of the 
highest luminosity sources. This improvement is 
responsible for most of the decrease in the mean 
Lboi from 5.3 L in Evans et al. (2009) to 4.3 L 
in this study. 

5. Our derived luminosity distribution is significantly 
different from that of Kryukova et al. (2012), who 
find a strong excess of sources at Lboi < 1-0 L 
compared to our results. Some of this discrepancy 
can be explained by the fact that we have modified 
the SEDs used to calculate Lboi, requiring a modi- 
fication to the Kryukova et al. relationship used to 
calculate L bo i, but some is also due to source selec- 
tion. By not requiring a (sub) millimeter detection 
they identify nearly 40% more protostars than we 
do in the same clouds, most of which are located 
at low Lboi- Future work is needed to better char- 
acterize the completeness of our sample and the 
reliability of the Kryukova et al. sample. 

6. 100 out of 230 protostars (43%) lack any available 
data in the far-infrared and submillimeter. The 
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calculated Lboi for these sources underestimates 
the true Lboi by a factor of 2.5 on average, and 
up to factors of 8 — 10 in extreme cases. Includ- 
ing far-infrared and submillimeter data for these 
sources once they become available from the Her- 
schel Gould Belt survey will likely increase both 
the mean and median Lboi by 35% - 40%. 

7. The conclusions of several recent theoretical studies 
of the protostellar accretion process that compare 
to the Evans et al. (2009) c2d luminosity distri- 
bution remain valid since our results are not sub- 
stantially different from the c2d results. As these 
studies demonstrate that there is more than one 
plausible accretion scenario that can match obser- 
vations, we have emphasized that future theoretical 
work is needed to better decipher the implications 
of protostellar luminosities for the underlying mass 
accretion process. 

We believe that the results presented here are the most 
complete and reliable census of protostellar luminosities 
assembled to date. Nevertheless, as outlined above in 
several avenues of future work must be pursued to 
better define the true completeness and reliability of the 
protostellar sample, more accurately measure Lboi for 
each source, and better understand the full extent and 
shape of the low end of the luminosity distribution. 
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TABLE 4 

List of Protostars and Basic Properties 



Spitzer 
Source Name 



Index 


Cloud 


(SSTc2d or SSTgb +) 


A v a 


1 


Aquila 


J1829053— 014157 


12.4 


2 


Aquila 


J1829234-013856 


12.4 


3 


Aquila 


J1829381— 015101 


12.4 


4 


Aquila 


J1829387-015100 


12.4 


5 


Aquila 


J1829419— 015012 


12.4 


6 


Aquila 


J1829433— 015652 


12.4 


7 


Aquila 


T1 OOn A 7H n 1 ET ET A O 

J 1829470— 015548 


12.4 


8 


Aquila 


J 1829595— 020106 


12.4 


9 


Aquila 


J 183U01 1 — 020609 


12.4 


10 


Aquila 


ti o o Ann c nonoco 

J 1830025—020258 


12.4 


11 


Aquila 


J 1830175— 020958 


12.4 


12 


Aquila 


JlooUz4o— U 15411 


12.4 


13 


Aquila 


JlooUzoy— (JzlU4o 


12.4 


14 


Aquila 


J 1830293— 015643 


12.4 


15 


Aquila 


J 1830469— 015646 


12.4 


16 


Aquila 


T1 Clin ft n 1 f-r*nri 

J 1830487— 015602 


12.4 


17 


Aquila 


Jloolozz— UzUlzo 


12.4 


18 


Aquila 


TI OOOI OO A1 FTOfl 

J 1832132— 015730 


12.4 


19 


Auriga/ California 


TA/I 1 A/I 1 A i OOAOnC 

J0410416+380805 


10.0 


20 


Auriga/California 


J0430036+351420 


7.5 


21 


Auriga/California 


TA /I OAAOA 1 'in IT A 

J0430082+351410 


7.5 


22 


Auriga/California 


ta>ioai ir i or'iooo 

J0430145+351332 


23.4 


23 


Cepheus 


J 2035463+675302 


5.4 


24 


Cepheus 


J2036198+675631 


6.5 


25 


Cepheus 


J2040567+672305 


5.4 


26 


Cepheus 


J2057130+773543 


5.4 


27 


Cepheus 


J2100207+681316 


5.4 


28 


Cepheus 


J2100221+681258 


5.4 


29 


Cepheus 


J2101328+681120 


5.4 


30 


Cepheus 


J2102212+675420 


5.4 


31 


Cepheus 


J2102273+675418 


5.4 


32 


Cepheus 


J2228030+690116 


5.4 


33 


Cepheus 


J2228074+690038 


5.4 


34 


Cepheus 


J2229333+751316 


5.4 


35 


Cepheus 


J2229594+751403 


5.4 


36 


Cepheus 


J2230318+751409 


5.4 


37 


Cepheus 


J2231056+751337 


5.1 


38 


Cepheus 


J2238428+751136 


5.4 


39 


Cepheus 


J2238469+751133 


5.4 


40 


Cepheus 


J2238530+751123 


5.4 


41 


Cepheus 


J2039062+680215 


5.4 


42 


Cepheus 


J2045539+675738 


5.4 


43 


Cepheus 


J2235234+751707 


5.4 


44 


Chamaeleon I 


J1104227- 771808 


6.6 


45 


Chamaclcon I 


Jl 106464— 772232 


6.6 


46 


Chamaeleon I 


J1106580-772248 


6.6 


47 


Chamaclcon I 


J1107161-772306 


6.6 


48 


Chamaclcon I 


J1107213— 772211 


6.6 


49 


Chamaclcon I 


J1107435— 773941 


5.9 


50 


Chamaeleon I 


J1108029-773842 


6.6 


51 


Chamaeleon I 


J1109285-763328 


6.6 


52 


Chamaclcon I 


J1109461-763446 


5.9 


53 


Chamaclcon I 


J1109472-772629 


9.2 


54 


Chamaeleon I 


J1110033-763311 


6.6 


55 


Chamaclcon I 


J1110113-763529 


8.0 


56 


Chamaclcon I 


J1111107— 764157 


6.6 


57 


Chamaeleon II 


J1253172-770710 


10.5 


58 


Chamaeleon II 


J1253428— 771511 


4.0 


59 


Chamaclcon II 


J1259065— 770739 


4.0 


60 


Corona Australis 


J1901480-365722 


7.9 


01 


Corona Australis 


J lyU14o4 — ODD f 14 


7 Q 


62 


Corona Australis 


J1901537-370033 


1.5 


63 


Corona Australis 


J1901585-365708 


7.9 


64 


Corona Australis 


J1902586-370735 


7.9 


65 


Corona Australis 


J1901086-365720 


7.9 


66 


Corona Australis 


J1903068-371249 


7.9 


67 


Corona Australis 


J1901506-365809 


7.9 


68 


Corona Australis 


J1901415-365831 


7.9 


69 


Corona Australis 


J1901553-365721 


7.9 


70 


Corona Australis 


J1901564-365728 


7.9 


71 


IC5146 


J2145585+473601 


3.6 


72 


IC5146 


J2147227+473214 


3.6 


73 


Lupus 


J1539277-344617 


1.0 



Observed Extinction Corrected 





Tbol 


-^bol 




Tbol' 


-i-bol' 




a 


(K) 


(Lq) 


a' 


(K) 


(L©) 


FIR/SI 


1.14 


80 


2.7 


0.96 


96 


2.9 


N 


1.39 


250 


0.84 


1.08 


320 


1.2 


K 


1.30 


240 


0.33 


1.09 


260 


0.46 


N 


0.82 


84 


3.0 


0.65 


110 


3.3 


N 


0.65 


110 


0.067 


0.46 


150 


0.078 


K 


0.81 


310 


0.67 


0.94 


380 


0.99 


K 


0.60 


220 


0.85 


0.29 


330 


1.1 


N 


0.74 


180 


0.037 


0.53 


230 


0.049 


N 


0.76 


110 


0.85 


0.34 


160 


0.98 


N 


0.48 


390 


1.2 


0.36 


470 


1.9 


N 


1.32 


130 


0.24 


1.11 


150 


0.30 


IN 


U. 21 


4UU 


n 7Q 


U.lz 


con 

oyu 


i.i 


V 


1 QQ 
l.OO 


yo 


O.O 


1 on 
1.2V 


11U 


a n 
4.U 


V 


U. 21 


4oU 


1 K 
1.0 


— U.UO 


71 n 
/ 1U 


z.O 


1\ 


1. lo 


i on 
12U 


U.10 


U.o 1 


14U 


n on 


V 


0.98 


190 


0.28 


0.80 


250 


0.36 


1\ 


1.00 


i fin 
10U 


U.lz 


1 Q K 

l.oO 


i 7n 
1 1 U 


U. 10 


IN 


n oq 


/ 


U.zo 


U. / o 


pro 

oy 


n oo 


1N 


— 0.32 


50 


7 n 
AO 


—0.88 


98 


7 Q 

A3 


■v 

I 


0.75 


98 


0.75 


0.80 


140 


0.80 


N 


0.46 


130 


0.81 


0.44 


180 


0.90 


N 


—0.39 


720 


0.94 


— 1.26 


1500 


4.7 


N 


1.00 


50 


1.4 


0.56 


53 


1.4 


Y 


—0.50 


760 


2.2 


—0.64 


2500 


4.9 


Y 


0.81 


150 


0.11 


0.30 


200 


0.12 


N 


0.10 


360 


2.0 


—0.04 


530 


2.4 


Y 


1.60 


89 


0.90 


1.36 


100 


0.95 


Y 


0.95 


87 


0.83 


0.52 


100 


0.86 


Y 


1.00 


21 


3.1 


0.94 


26 


3.1 


Y 


1.10 


180 


0.41 


0.64 


210 


0.44 


Y 


—0.29 


36 


0.32 


—0.64 


47 


0.33 


Y 


0.84 


160 


2.1 


0.40 


180 


2.3 


Y 


1.00 


42 


0.94 


0.90 


43 


0.95 


Y 


0.12 


210 


0.086 


0.10 


240 


0.096 


N 


0.49 


270 


0.32 


0.16 


330 


0.36 


N 


0.71 


35 


0.43 


0.64 


36 


0.43 


Y 


—0.51 


30 


0.30 


—0.77 


31 


0.30 


Y 


1.20 


72 


1.8 


0.89 


78 


1.8 


Y 


1.60 


110 


6.2 


1.94 


120 


6.6 


Y 


0.97 


90 


1.4 


0.44 


150 


1.5 


Y 


0.94 


35 


4.0 


0.79 


35 


4.1 


Y 




340 


48 




450 


62 


Y 


2.40 


230 


8.7 


3.69 


280 


9.5 


Y 


0.61 


350 


0.086 


0.27 


510 


0.11 


N 


1.23 


66 


0.69 


1.11 


74 


0.72 


Y 


0.51 


170 


0.15 


0.92 


210 


0.17 


N 


—0.10 


490 


0.076 


—0.03 


590 


0.099 


N 


—0.08 


650 


0.089 


—0.45 


820 


0.13 


N 


— 1.05 


1400 


0.45 


—1.43 


1700 


0.85 


N 


—0.10 


710 


0.97 


—0.39 


900 


1.4 


N 


1.17 


260 


0.90 


1.30 


300 


1.1 


Y 


—0.41 


720 


0.18 


—1.13 


1100 


0.26 


N 


—0.81 


1100 


0.092 


—1.25 


1500 


0.21 


N 


0.32 


270 


0.018 


—0.18 


430 


0.022 


N 


—0.44 


1100 


0.37 


— 1.24 


1600 


0.80 


N 


0.09 


330 


0.0076 


—0.20 


650 


0.0095 


N 


—0.72 


660 


30 


—0.63 


1500 


63 


N 


0.65 


130 


0.43 


0.32 


160 


0.45 


Y 


0.68 


230 


1.7 


1.12 


260 


1.8 


Y 


0.78 


93 


4.4 


1.03 


130 


4.7 


Y 


1.41 


17 


3.7 


1.20 


23 


3.7 


Y 


— 1.09 


19 


1.4 


—1.29 


23 


1.4 


Y 


0.88 


13 


6.9 


0.72 


15 


7.0 


Y 


1.66 


61 


1.2 


1.48 


66 


1.3 


Y 


-0.80 


1000 


3.5 


-1.22 


1500 


6.9 


Y 


0.36 


460 


8.8 


-0.21 


610 


13 


Y 


0.92 


210 


15 


1.01 


270 


17 


Y 


0.75 


270 


7.1 


0.27 


390 


8.5 


Y 


2.64 


260 


0.44 


2.47 


310 


0.53 


Y 


2.78 


200 


1.4 


2.56 


210 


1.8 


Y 


0.82 


150 


7.4 


1.11 


170 


7.8 


Y 


0.74 


86 


35 


0.63 


90 


36 


Y 


-0.84 


2700 


0.94 


-1.20 


3300 


1.3 


N 
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Spitzer 






Observed 




Extinction Corrected 








Source Name 






bol 


r 

-^bol 




Tbol' 


ibol' 


FTR /SMM b 


T H 

n cx 


Cloud 


(Jul e^u or ooigu ^ j 








l^Q ; 


a' 


(K) 


(L©) 


74 


Lupus 


J1539282-344618 


2.0 


— 0.84 


2100 


0.24 


— 1.11 


2600 


0.36 


N 


75 


Lupus 


J1607100-391103 


3.0 


— 1 .04 


2100 


0.69 


— 1.23 


3700 


1.6 


N 


76 


Lupus 


J1608217-390421 


1.0 


-1.13 


2600 


0.21 


—1.46 


3100 


0.28 


N 


77 


Lupus 


J1608224-390446 


0.0 


-0.48 


1900 


1.7 


—0.63 


1900 


1.7 


N 


78 


Lupus 


J1609180-390453 


2.9 


1.10 


39 


0.41 


1.14 


39 


0.41 


Y 


79 


Ophiuchus 


J1625381-242236 


13.2 


-0.79 


1100 


0.11 


— 1.22 


1600 


0.36 


N 


80 


Ophiuchus 


J1625561-242048 


4.5 


-0.60 


1100 


0.89 


—0.88 


1400 


1.3 


N 


81 


Ophiuchus 


J1626103-242054 


27.0 


-0.46 


290 


1.0 


— 1.43 


1300 


3.7 


Y 


82 


Ophiuchus 


J1626146-242507 


0.0 




7 


0.034 




7 


0.034 


N 


83 


Ophiuchus 


J1626188-242819 


19.7 


-0.73 


990 


0.50 


—1.28 


1600 


2.5 


N 


84 


Ophiuchus 


J1626213-242304 


9.8 


1.46 


210 


8.6 


1.14 


250 


11 


Y 


85 


Ophiuchus 


J1626236-244314 


4.0 


-1.12 


1500 


0.48 


— 1.34 


1700 


0.73 


Y 


86 


Ophiuchus 


J1626240-241613 


13.3 


-0.71 


980 


1.9 


— 1.09 


1500 


5.5 


N 


87 


Ophiuchus 


J1626254- 242301 


9.8 


0.87 


140 


0.010 


0.60 


200 


0.12 


N 


88 


Ophiuchus 


J1626256-242428 


9.8 


1.65 


72 


0.038 


1.44 


84 


0.043 


N 


89 


Ophiuchus 


J1626404- 242714 


9.8 


0.45 


380 


0.064 


0.32 


470 


0.092 


N 


90 


Ophiuchus 


J1626441-243448 


9.8 


2.49 


330 


0.98 


2.50 


380 


1.4 


N 


91 


Ophiuchus 


J1626450-242307 


18.5 


-0.64 


820 


0.29 


— 1.20 


1500 


1.2 


N 


92 


Ophiuchus 


J1626484- 242838 


9.8 


0.02 


440 


0.12 


—0.06 


570 


0.17 


N 


93 


Ophiuchus 


J1626584- 244531 


10.7 


-0.45 


840 


1.3 


—0.70 


1300 


2.7 


N 


94 


Ophiuchus 


J1627023-243727 


9.8 


1.53 


420 


3.3 


0.91 


520 


4.8 


N 


95 


Ophiuchus 


J1627029-242614 


9.8 


-0.19 


380 


0.036 


—0.23 


550 


0.050 


N 


96 


Ophiuchus 


J1627052-243629 


9.8 


1.27 


97 


0.16 


1.07 


120 


0.18 


N 


97 


Ophiuchus 


J1627067- 243814 


9.8 


0.61 


330 


0.48 


0.73 


420 


0.64 


N 


98 


Ophiuchus 


J1627094- 243718 


9.8 


1.69 


370 


13 


1.06 


420 


18 


N 


99 


Ophiuchus 


J1627158-243843 


18.6 


-0.70 


220 


0.63 


— 1.64 


830 


1.1 


N 


100 


Ophiuchus 


J1627175-242856 


9.8 


0.25 


190 


0.52 


0.37 


260 


0.61 


N 


101 


Ophiuchus 


J1627214- 244143 


9.8 


-0.03 


610 


1.1 


—0.02 


720 


1.8 


N 


102 


Ophiuchus 


J1627218-242727 


9.8 


-0.05 


180 


0.019 


—0.08 


290 


0.022 


N 


103 


Ophiuchus 


J1627245-244103 


9.8 


1.01 


170 


0.31 


1.23 


230 


0.37 


N 


104 


Ophiuchus 


J1627269-244050 


9.8 


1.17 


240 


2.6 


1.15 


300 


3.3 


N 


105 


Ophiuchus 


J1627279-243933 


9.8 


2.29 


260 


5.0 


2.13 


280 


7.1 


N 


106 


Ophiuchus 


J1627284- 242721 


9.8 


-0.03 


310 


0.48 


—0.13 


450 


0.63 


N 


107 


Ophiuchus 


J1627301-242743 


9.8 


-0.12 


500 


0.97 


—0.02 


620 


1.5 


N 


108 


Ophiuchus 


J1627372-244237 


9.8 


0.13 


460 


0.12 


0.15 


560 


0.18 


N 


109 


Ophiuchus 


J1627398-244315 


9.8 


-0.15 


570 


0.72 


—0.15 


690 


1.1 


N 


110 


Ophiuchus 


J1628216-243623 


9.8 


1.23 


33 


0.24 


0.96 


36 


0.24 


Y 


111 


Ophiuchus 


J1628578-244054 


9.8 


0.67 


320 


0.027 


0.58 


430 


0.037 


N 


112 


Ophiuchus 


J1631356-240129 


9.8 


0.14 


270 


1.6 


0.14 


390 


2.0 


Y 


113 


Ophiuchus 


J1631367-240419 


9.8 


-0.27 


74 


0.17 


—0.22 


160 


0.19 


Y 


114 


Ophiuchus 


J1631437-245524 


9.8 


0.23 


520 


0.26 


0.19 


690 


0.40 


N 


115 


Ophiuchus 


J1631520-245726 


9.8 


0.82 


120 


0.0082 


0.61 


150 


0.0095 


N 


116 


Ophiuchus 


J1631524- 245536 


9.8 


1.07 


260 


0.11 


0.87 


330 


0.15 


N 


117 


Ophiuchus 


J1632009-245642 


9.8 


1.39 


140 


2.5 


1.36 


180 


2.8 


Y 


118 


Ophiuchus 


J1632226-242831 


9.8 


5.03 


45 


8.5 


4.87 


45 


8.8 


Y 


119 


Ophiuchus 


J1633556-244205 


4.6 


-1.22 


1500 


0.17 


—1.31 


1800 


0.29 


N 


120 


Ophiuchus North 


J1646582-093519 


5.6 


0.66 


230 


0.44 


0.61 


280 


0.50 


Y 


121 


Ophiuchus North 


J1648456- 141636 


5.6 


-0.97 


1400 


1.2 


— 1.39 


1700 


2.2 


N 


122 


Ophiuchus North 


J1657196- 160923 


5.7 


2.40 


38 


0.80 


2.57 


39 


0.82 


Y 


123 


Perseus 


J0325223+304513 


5.9 


2.34 


52 


2.0 


2.20 


53 


2.1 


Y 


124 


Perseus 


J0325362+304515 


5.9 


1.59 


12 


8.5 


1.49 


12 


8.5 


Y 


125 


Perseus 


J0325364+304522 


5.9 


2.62 


66 


4.9 


2.55 


70 


5.2 


Y 


126 


Perseus 


J0325388+304406 


5.9 


2.16 


47 


6.9 


2.04 


48 


7.0 


Y 


127 


Perseus 


J0325391+304358 


5.9 


2.36 


160 


0.69 


2.02 


170 


0.78 


N 


128 


Perseus 


J0326374+301528 


5.9 


1.09 


64 


0.91 


1.02 


73 


0.94 


Y 


129 


Perseus 


J0327382+301358 


5.9 


-0.19 


260 


0.78 


—0.40 


360 


0.89 


Y 


130 


Perseus 


J0327390+301303 


5.9 


2.68 


62 


3.5 


2.45 


65 


3.6 


Y 


131 


Perseus 


J0327432+301228 


5.9 


2.39 


54 


1.7 


2.23 


( 


1 *7 
1. / 


Y 


132 


Perseus 


J0327476+301204 


5.9 


-0.09 


740 


2.5 


—0.25 


950 


3.6 


Y 


133 


Perseus 


J0328003+300801 


5.9 


0.95 


230 


0.25 


0.99 


280 


0.29 


Y 


134 


Perseus 


J0328325+311105 


5.9 


0.78 


52 


0.26 


0.50 


74 


0.27 


Y 


135 


Perseus 


J0328344+310051 


5.9 


0.83 


240 


1.1 


0.88 


290 


1.3 


Y 


136 


Perseus 


J0328345+310705 


5.9 


0.54 


150 


0.12 


0.32 


170 


0.13 


Y 


137 


Perseus 


J0328350+302009 


5.9 


0.15 


49 


0.36 


—0.22 


89 


0.36 


Y 


138 


Perseus 


J0328370+311330 


5.9 


2.35 


100 


9.5 


1.94 


110 


10 


Y 


139 


Perseus 


J0328391+310601 


5.9 


1.68 


28 


0.23 


1.55 


29 


0.23 


Y 


140 


Perseus 


J0328397+311731 


5.9 


0.57 


250 


0.18 


0.61 


300 


0.21 


Y 


141 


Perseus 


J0328406+311756 


5.9 


1.02 


12 


0.58 


0.95 


12 


0.58 


Y 


142 


Perseus 


J0328432+311732 


5.9 


0.36 


490 


1.6 


0.25 


640 


2.0 


Y 


143 


Perseus 


J0328453+310541 


5.9 


1.11 


62 


0.41 


1.08 


72 


0.43 


Y 


144 


Perseus 


J0328555+311436 


5.9 


3.03 


54 


22 


2.37 


55 


22 


Y 


145 


Perseus 


J0328563+312227 


5.9 


-0.14 


440 


0.15 


-0.40 


620 


0.19 


N 


146 


Perseus 


J0328573+311415 


5.9 


1.60 


100 


5.3 


1.40 


110 


5.7 


Y 


147 


Perseus 


J0328584+312217 


5.9 


0.83 


240 


0.96 


0.83 


280 


1.1 


Y 
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Spitzer Observed Extinction Corrected 

Source Name T bol L bol T bol ' L bol ' 



Index 


Cloud 


(SSTc2d or SSTgb +) 


A v a 


a 


(K) 


(L©) 


a.' 


(K) 


(L©) 


FIR/SMM b 


1 A Q 

14o 


Perseus 


JUozo0yo-|-oll04o 


n 

o.y 


n nft 
U.Uo 


1 n 
1U 


Q O 
O.Z 


n no 
— U.Uo 


1 /in 
14U 


Q A 
0.4 


V" 

1 


149 


Perseus 


TnQonnnc; i qi 1 onn 
JUoZyUUo+ollzUU 


5.9 


2. 16 


30 


0.66 


1.96 


30 


0.67 


Y 


150 


Perseus 


Tnuonni c; i qi onon 
JUoZyUlO+olzUZU 


5.9 


2.09 


230 


8.2 


2.30 


270 


9.4 


■v 
Y 


151 


Perseus 


TnQonnQQ i qi oqi a 
JUoZyUoo+olZol4 


5.9 


1.13 


320 


0.088 


1.01 


370 


0.11 


1N 


152 


Perseus 


TnQonnQ'7 i qi 1 en 1 ! 
JUozyUo / +0II0U0 


5.9 


1.21 


170 


33 


0.85 


220 


07 
0/ 


Y 


153 


Perseus 


TnQonn/in 1 q 1 1 A A a 
J UozyU4U+o 1 1440 


n 

o.y 


1.43 


1 1 
1 / 


0.68 


1.31 


19 


0.68 


■v 
Y 


1 C A 

104 


Perseus 


TflQOnn77 1 QI 01 £."7 

juozyu/ /+01210/ 


c: n 

o.y 


2. 18 


OQn 
zoU 


18 


2.31 


oftn 
zoU 


01 

Zl 


Y 


loo 


Perseus 


TnQom n/i 1 qi 1 qqi 
JUozyiU4+ollool 


c; n 

o.y 


z.Oo 


Q "I 

ol 


1 n 

/ .y 


O A Q 

z.4o 


Q 1 
Ol 


q n 
o.U 


Y 


156 


Perseus 


TAQoni nft 1 qi 1 con 
JUozyiUo+olloZU 


5.9 


1.95 


56 


3.2 


1.84 


58 


3.3 


Y 


157 


Perseus 


TnQomio 1 QI 1 OQI 
JUoZyilZ-T-ollool 


5.9 


1.94 


29 


1.2 


1.77 


32 


1.3 


■V" 

1 


1 e.q 
IOo 


Perseus 


TnQom on 1 qi 1 qhc; 
JUoZyiZU-holloUO 


n 

o.y 


n c\q 

u.yo 


zO 


A O 
4.Z 


n q 1 
U.ol 


oc; 
ZO 


A O 
4.Z 


Y 


159 


Perseus 


TnQom on i qi 1 01 a 
JUoZyiZy+ollol4 


5.9 


1.05 


240 


0.100 


1.58 


270 


1.1 


1N 


160 


Perseus 


TnQom q c; i qi 1 oco 
JUoiyioO+olloOO 


5.9 


2.41 


35 


0.84 


2.50 


36 


0.85 


Y 


101 


Perseus 


T0Q00171 I QlOV/lft 

juozyi / 1+0 iz i 40 


k n 

o.y 


1.(0 


QO 
OZ 


n ftp; 
U.OO 


1 7n 
1. / U 


QA 

o4 


n ftft 
U.OO 


Y 


1 ftO 

10Z 


Perseus 


TnQom co 1 QioQin 

jUv5zy loz+oizoiy 


k n 

o.y 


1 Oft 

l.zo 


Ol 

zl 


U.04 


111 
1.11 


OA 

Z4 


n K/i 
U.04 


V" 

I 


163 


Perseus 


THQOni C7 1 QIOQOCT 

JUoiyio / -\-6YZ6Zb 


5.9 


—0.23 


190 


1.9 


—0.67 


410 


2.1 


\7 
I 


164 


Perseus 


TnQononn 1 qi O/inv 
JUozyzUU+ol24U / 


5.9 


0.42 


75 


1.5 


0.34 


10 


1.6 


■V" 
I 


165 


Perseus 


TnQOnOQ/1 1 Q1 QQOn 

JUozyZo4-|-oloozy 


5.9 


1.51 


60 


0.36 


1.51 


ft A 

o4 


0.37 


■V" 
I 


166 


Perseus 


TnQonci 1 qi Qnnft 
JUozyolo+oloyUo 


c n 

o.y 


Q A A 
0.44 


q n 
oy 


n c;n 

U.oU 


Q AC\ 

o.4U 


a n 
4U 


0.51 


\7 
I 


1 ft7 

lo 7 


Perseus 


TnQQm £i 1 QnoQ/in 
JUooUlol+oUZv54y 


cc n 

o.y 


1.70 


93 


1 A 

1.4 


1 CC7 
1.0 / 


100 


1.5 


Y 


1 ftQ 

IOo 


Perseus 


TnQQnQOA 1 QnoftOft 
JUooUoZ0-|-oUZ0Z0 


c; n 

o.y 


no 
z.Uo 


QA 


n 1 ft 
U.10 


70 
z. / z 


Q 

oO 


n 1 ft 
U. 10 


V" 

1 


169 


Perseus 


TnQQ 1 onn i on a r; on 

JUoolzUy+oU4ooU 


5.9 


0.98 


32 


1.2 


1.48 


32 


1.2 


\j 
1 


170 


Perseus 


TnQQOI 7 7 n 1 QC\AC\A'7 

JUoozl /y-|-oU4y4/ 


5.9 


1.07 


25 


1.3 


0.87 


26 


1.3 


■v 
Y 


i/i 


Perseus 


TnQQOom 1 Qino/in 
J UooZZy 1+01UZ4U 


x. n 

o.y 


n An 
U.4U 


1 on 
IzU 


n c;o 
U.Oz 


n 1 ft 
U.10 


1 ftn 
10U 


n c^ft 
U.OO 


V" 

1 


1 70 
1 / Z 


Perseus 


TnQQQnoc; i Qinc".Qi 
JUoooUy0-+-olU0ol 


n 

o.y 


1 1 Q 

1. lo 


01 n 
Z1U 


n c\a 
U.U4z 


n ftn 

u.oy 


Qon 
oZU 


n c\a n 
U.U4y 


TvT 

IN 


173 


Perseus 


TnQQQIOQ 1 Q1010/1 

JUooolZo-|-olZlZ4 


5.9 


0.41 


480 


2.9 


0.16 


610 


3.8 


Y 


174 


Perseus 


TnQQQi qo i qi onnc; 
JUoooloo+olZUUO 


5.9 


1.44 


59 


0.100 


1.31 


69 


0.11 


Y 


1 

1 /O 


Perseus 


TnQQQi/iQ 1 Qiri7in 
JUoool4o+olU / 1U 


5.9 


2.22 


34 


0.63 


2.73 


36 


0.64 


Y 


176 


Perseus 


TnQQQift/i 1 Qmcco 
JUooolo44-Jl(Jooz 


5.9 


1.73 


26 


1.0 


1.56 


26 


1.0 


■v 
Y 


177 


Perseus 


JU3o31oo+olU755 


5.9 


1.57 


100 


1.6 


1.69 


120 


1.7 


Y 


1 7Q 
I/O 


Perseus 


TOQQQ1 7Q 1 QinnQI 

JUoool /o+olUyol 


K n 

o.y 


Q QQ 
0.00 


A 1 
41 


Q c; 
0.0 


no 

z.yz 


A A 

44 


Q ft 
O.O 


V" 

1 


1 TO 

i / y 


Perseus 


TnQQQOnQ I Qiri701 

J{Jooo2Uo-\-o 1U / Zl 


t; n 

o.y 


n qq 
U.00 


A 1 
4 / 


n ftn 
U.OU 


n fti 
U.01 


00 


n fti 
U.01 


V" 

I 


180 


Perseus 


TnQQQOvo i Qin7in 
JU000Z / Z+olU / 1U 


ec n 

o.y 


1.93 


ftO 
OZ 


1.3 


1.79 


ftft 
OO 


1 A 

1.4 


Y 


1 Q 1 

lol 


Perseus 


TnQ/ionoi 1 Qi/iQno 
JUo4ZUZl-(-ol4oUZ 


t; n 

o.y 


1 A 1 
1.4 / 


i nn 

iyu 


n nftn 

u.uoy 


1 OQ 

l.zo 


o/in 
Z4U 


n non 
U.UoU 


IN 


loZ 


Perseus 


TnQ/iQ/ic".i 1 QonQt;Q 
J Uo4o40 1 -hozUoOo 


x. n 

o.y 


n 00 
— U.zz 


04U 


n 71 
U. / 1 


n 1 
U.lz 


ft c^n 
OOU 


n oq 

u.yo 


IV 
1N 


183 


Perseus 


TnQ a oenn 1 QonQO a 
JUo4ooUy+ozUoz4 


5.9 


1.51 


16 


0.68 


1.38 


17 


0.69 


Y 


184 


Perseus 


TnQ/iQp;in i QonQno 
JUo4oolU+ozUoUo 


5.9 


—0.28 


52 


0.40 


—0.30 


57 


0.41 


Y 


1 Q f"» 
loO 


Perseus 


TnQ a Q^K^j-Qonn^o 
JUo4o000-t-oZUU0Z 


o.y 


n c;n 
U.OU 


OQ 

zo 


1 ft 
1.0 


n 

U.OO 


OQ 

Zo 


1 ft 
1.0 


v 

1 


1 Ofi 
loO 


Perseus 


inQ/iQtrfio 1 QonQn/i 
JUo4o00o-t-oZUoU4 


K n 

o.y 


1 Q*7 

1.0 / 


00 

zz 


1 A 
1.4 


1 QQ 
l.OO 


OQ 

Zo 


1.4 


Y 


187 


Perseus 


TnQ AOKaa. 1 Qom ka 
JUo4ooyo+ozUlo4 


19.8 


—0.34 


620 


1.8 


—0.90 


1300 


6.2 


V" 

1 


1 QQ 

loo 


Perseus 


TnQ/i/ino/i 1 Qonon/i 
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-0.08 


700 


64 


Y 


221 


Serpens 


J1829578+011251 


9.6 
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TABLE 4 — Continued 







Spitzer 






Observed 




Extinction Corrected 








Source Name 
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0.46 
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N 
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0.18 
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N 



a Value of Ay used for dereddening, as explained in the text. 

k Flag indicating that the protostar docs ( Y) or docs not (N) have at least one observed photometry point in the far-infrared or submillimeter 
(70 fim < A < 850 /mi; sec £|4~3l for details). 



